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Autism spectrum disorder (ASD) is a complex heterogeneous neurodevelopmental 
disorder characterized by alterations in social functioning, communicative abilities, and 
engagement in repetitive or restrictive behaviors. The process of aging in individuals with 
autism and related neurodevelopmental disorders is not well understood, despite the fact 
that the number of individuals with ASD aged 65 and older is projected to increase by 
over half a million individuals in the next 20 years. To elucidate the effects of aging in 
the context of a modified central nervous system, we investigated the effects of age 
on the BTBR T + tf/j mouse, a well characterized and widely used mouse model that 
displays an ASD-like phenotype. We found that a reduction in social behavior persists into 
old age in male BTBR T + tf/j mice. We employed quantitative proteomics to discover 
potential alterations in signaling systems that could regulate aging in the BTBR mice. 
Unbiased proteomic analysis of hippocampal and cortical tissue of BTBR mice compared 
to age-matched wild-type controls revealed a significant decrease in brain derived 
neurotrophic factor and significant increases in multiple synaptic markers (spinophilin, 
Synapsin I, PSD 95, NeuN), as well as distinct changes in functional pathways related 
to these proteins, including " Neural synaptic plasticity regulation" and "Neurotransmitter 
secretion regulation." Taken together, these results contribute to our understanding of 
the effects of aging on an ASD-like mouse model in regards to both behavior and protein 
alterations, though additional studies are needed to fully understand the complex interplay 
underlying aging in mouse models displaying an ASD-like phenotype. 
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INTRODUCTION 

Autism spectrum disorder (ASD) is a heterogeneous neu- 
rodevelopmental disorder characterized by varying degrees of 
altered social functioning and engagement in repetitive, stereo- 
typed behaviors accompanied by display of narrowed inter- 
ests (American Psychiatric Association, DSM-V 2012). The 
most recent surveillance data collected in 2008 by the Center 
for Disease Control's Autism and Developmental Disabilities 
Monitoring Network indicates that approximately 1 in 88 chil- 
dren has an autism spectrum disorder (CDC, 2010). The United 
States Census Bureau in 2006 projected a doubling of the US pop- 
ulation aged 65 and older by 2030. Assuming the life expectancy 
of individuals with ASD is similar to that of the general popula- 
tion, based on current prevalence rates in school aged children, 
this population increase would result in a prevalence of approxi- 
mately 700,000 individuals with ASD aged 65 in the next 20 years 
(Piven and Rabins, 2011). 



Abbreviations: BDNF, brain derived neurotrophic factor; BTBR, black and tan, 
brachyuric; ASD, autism spectrum disorder. 



Despite the population statistics indicating a significant num- 
ber of aged individuals with ASD, studies that investigate the 
effects of the neurodevelopmental disabilities including ASD, in 
an aged physiological context are lacking. Due to the paucity and 
limitations of studies and the expected increase in the aged indi- 
viduals presenting ASD a multidisciplinary conference was held 
in March 2010 to identify gaps in knowledge regarding older 
adults with ASD. The panel concluded that descriptive stud- 
ies of aging in humans and rodent models were needed (Piven 
and Rabins, 2011). A limited number of human studies have 
investigated the persistence of the ASD core deficits into adult- 
hood (Ballaban-Gil et al., 1996; Seltzer et al., 2003; Shattuck 
et al., 2007), childhood factors that are associated with prognosis 
(Rumsey et al., 1985), and outcome studies (Rutter et al., 1967; 
Howlin et al, 2004; Renty and Roeyers, 2006). Major findings 
from these studies suggest that language development (Kanner 
et al., 1972) and Intelligence Quotient (IQ) predicts outcome 
(Lotter, 1974; Rumsey et al, 1985; Szatmari et al, 1989). The 
ASD diagnosis persists with inconsistent findings about which, 
if any, core deficits change (Ballaban-Gil et al., 1996; Piven et al., 
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1996; Seltzer et al, 2003), and that the adults with ASD remain 
dependent on others (Rutter et al, 1967) even if they have a 
normal IQ (Howlin et al, 2004). Despite persistent dependence 
some studies found that individuals with ASD had improved 
symptoms over time (Szatmari et al., 1995; Piven et al, 1996; 
Seltzer et al., 2003). However, a few limitations hinder these 
early studies of aging in ASD, including study individuals being 
only 20-30 years old, small sample sizes, and an inability to 
completely generalize the study findings given the heterogene- 
ity of the population examined (Rumsey et al., 1985; Gillberg 
and Wing, 1999; Fombonne et al., 2001; Howlin and Moss, 
2012). 

Given the heterogeneous nature of ASD, mouse models have 
proven to be particularly useful and reliable in elucidating the 
mechanism of these disorders (Crawley, 2007; Moy and Nadler, 
2008). There are several approaches to ASD mouse models: 
targeted gene mutations, defects in neurotransmitters or brain 
regions, inbred strains, and models based on the comorbidity of 
other human diseases with autism (Crawley, 2007). One strain in 
particular, BTBR T + tf/j, has been shown to be an especially rele- 
vant animal model of ASD (Moy et al., 2007). Extensive behav- 
ioral characterization of the BTBR (black and tan, brachyuric) 
mouse model has revealed low sociability compared to wild- 
type strains (Bolivar et al, 2007; Moy et al., 2007), resistance 
to change (Moy et al., 2007; Moy and Nadler, 2008), increased 
display of self-grooming behavior (Pobbe et al., 2010), display 
of repetitive behaviors (Pearson et al., 2011), and reduced dis- 
play of territorial scent marking (Wohr et al, 201 1). Furthermore, 
unusual vocalizations have also been extensively characterized in 
BTBR mice (Scattoni et al., 2008, 2011), in addition to instances 
of social avoidance and gaze aversion (Defensor et al, 2011). 
While the autistic-like behavioral phenotype of the young male 
BTBR mouse has been studied extensively, the effects of age on 
this ASD-like mouse model have not been elucidated. As aged 
human autistic tissue samples are limited and the BTBR mouse 
model is well characterized and widely accepted, we selected this 
mouse model to study the nature of the ASD phenotype in an 
aged context. To enhance our understanding of the aged BTBR 
mouse and to gain further insight into the effects of aging of 
ASD, we followed our extensive behavioral characterization of 
the aged BTBR mouse with quantitative proteomic analyses on 
cortical and hippocampal tissues as these two brain regions have 
been strongly associated with ASD (Mundy, 2003; Nadler et al., 
2006). We hypothesized that both the behavioral phenotype and 
proteomic profile would be different between the aged ASD-like 
mice [or mice with an abnormal central nervous system (CNS)] 
and the aged wildtype mice (or mice with a presumed normal 
CNS). Our results contribute to our understanding of the effects 
of aging on an ASD-like mouse model, and it is our hope that 
further research will enhance our understanding of aging in an 
abnormal CNS. 

MATERIALS AND METHODS 
EXPERIMENTAL ANIMALS 

All experimental animal procedures were approved by the Animal 
Care and Use Committee of the National Institute on Aging, 
National Institutes of Health. Male BTBR T + tf/J mice (n = 8, 15 



months of age) and male control wildtype ( WT) C57BL6J (« = 5, 
15 months of age) were housed in the NIA animal facility on 
a 12 h light and dark cycle from 6 am to 6 pm and received ad 
libitum access to food and water throughout the duration of the 
study. All behavioral testing described in detail below was per- 
formed between the hours of 9 AM and 5 PM. Prior to behavioral 
testing, animals were moved from the vivarium to an experimen- 
tal room and given a minimum of 30 min to habituate to the new 
environment. Between each experimental animal, all apparatuses 
were wiped down with 70% ethanol followed by D. I. water to pre- 
vent any scent retention that could possibly bias the performances 
of later mice. 

SOCIAL PREFERENCE ASSESSMENT 

Social behavior was assessed using a three-chambered sociability 
apparatus with a protocol modeled after previously established 
methods (Crawley, 2007). In brief, animals were first habituated 
for 10 min to the three-chambered sociability apparatus (60 x 
40 x 22 cm, Stoelting, Wood Dale, IL) free of any stimuli and 
returned to their home cage for a minimum of 30 min. Prior to 
the return of the experimental mouse to the chamber for testing, 
another mouse of the same strain, age, and gender and an inani- 
mate object both housed in identical wire mesh cylinders (15 cm 
tall, 7 cm in diameter) were placed in opposing compartments of 
the apparatus. Experimental animals were then returned to the 
center chamber of the apparatus and were free to choose to inter- 
act with either the mouse or object enclosed in the respective wire 
mesh containers. Interactions were recorded for 10 min and ana- 
lyzed with ANY-Maze Video Tracking Software (Stoelting). A dis- 
crimination index (DI) was calculated that takes into account the 
relative exploration time spent with either the contained mouse or 
object: (animal time/total time — object time/total time) x 100. 
As socially-normal mice show a robust preference for another 
mouse over an inanimate object, a high DI is reflective of socially- 
normal behavior while a low DI is indicative of deficits in social 
functioning. 

NOVEL OBJECT PREFERENCE TEST 

The novel object preference test was performed as described pre- 
viously (Dere et al., 2007; Park et al., 2010). In brief, mice were 
given 10 min to explore a 25 cm 2 opaque walled box with two 
identical objects placed equidistant from each other in the cen- 
ter of the box. After 10 min of habituation to the two identical 
objects was completed, animals were returned to their home 
cages for a minimum of 30 min. After the habituation task, 
trained experimenters replaced one of the identical objects with 
a novel object that had never been seen by the experimental 
mouse. Experimental animals were then returned to the opaque 
box and given 10 min in which they were free to interact with 
either the familiar or novel object. Both the habituation session 
and testing session were videotaped and the animal's behavior 
automatically tracked by ANY-Maze Video Tracking Software 
(Stoelting). Discrimination indices were again calculated: (novel 
time/total time-same time/total time) x 100. A high discrimi- 
nation index (DI) is indicative of an enhanced ability to pref- 
erentially recognize the novel object, suggesting intact memory 
functioning. 
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FIGURE 1 | Social function, cognitive ability, general behavior, anxiety 
and motor coordination in an aged ASD mouse model. Social preference 
was assessed using a three-chambered sociability apparatus in aged male 
BTBR T + tf/j mice and age-sex-matched control mice (A, B). Percentage of 
time spent exploring both the contained inanimate object and the contained 
mouse was measured for male BTBR (n = 8) and control mice (n = 5) (A). A 
discrimination index (Dl) was calculated (B) which represents the degree to 
which experimental or control animals displayed a preference for exploring 
the contained mouse over the contained inanimate object. Cognitive 
functioning was assessed via the novel object preference task (NOP) (C, D). 
Percentage of time spent exploring either a familiar or a novel object was 
again measured for male BTBR (n = 8) and control mice (n = 5) (C) and a Dl 
was calculated (D) with a higher Dl indicating an increased ability to discern a 
novel object from a familiar object. Exploration of an open field apparatus was 



assessed in male BTBR and control mice was measured with regard to the 
following parameters: time spent exploring the peripheral or central areas of 
the apparatus relative to the total amount of time spent exploring (E), total 
distance traveled (F), number of ambulatory episodes (G), and number of 
jumps (H). Exploration of the open field apparatus was repeated with a 
Light/Dark insert add-on and relative time spent in either the light half or the 
dark half of the chamber was recorded (I). General anxiety was assessed 
using the elevated plus maze task. For both BTBR T + tf/j mice and age-and 
sex-matched controls, relative spent in the closed vs. the open arm was 
recorded as a percentage of total time (J). Finally, motor coordination was 
assessed using a rotarod apparatus latency to fall was recorded in both BTBR 
T + tf/j and control mice (K, n = 5). Data are expressed as means ± s.e.m. 
Asterisks represent p-values as shown: *p < 0.05, **p < 0.01, ***p < 0.001. 
Statistical significance was measured using a Student's t-test. 



OPEN FIELD TEST 

A generalized view of the behavioral characteristics of both 
experimental and control mice was obtained via an open field 
test performed as described previously (Holmes et al., 2002). 
In brief, each mouse was placed in a square chamber approx- 
imately 0.14 m 2 in size (Med Associates Inc, St. Albans, VT) 
and given 10 min to freely explore. Activity Monitor software 
(Med Associates Inc) automatically recorded multiple parame- 
ters over the duration of the 10 min test. Individual data points 
were then compiled and subsequently analyzed to produce the 
following measurements: distance traveled, time spent in the 
peripheral areas of the chamber, time spent in the center area 



of the chamber, number of ambulatory episodes, and number of 
jumping episodes. Relative time spent in the periphery vs. the cen- 
ter was calculated for each mouse by dividing the time spent in 
either the periphery or center by the total exploration time (600 s 
for all animals). 

LIGHT/DARK EXPLORATION TEST 

The light-dark exploration test was performed as previously 
described (Chadwick et al, 201 1) as a means to assess general anx- 
iety behavior. The same square chamber used in the Open field 
test was again used to complete the Light/Dark exploration test, 
with the additional insertion of a dark chamber taking up 50% 
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FIGURE 2 | Protein expression in wild type and BTBR mouse 
cortex. Western blotting was performed using mouse cortical protein 
lysates (A). Protein levels were measured in BTBR and control mice 
with respect to the following proteins in the cortex: pro-BDNF (B), 
mature BDNF (C), phospho-TrkB (D), total TrkB (E), Akt (F) 



phospho-synapsin 1 (G), synapsin 1 (H) synaptophysin (I), PSD95 (J), 
spinophilin (K), Neuronal cell marker NeuN (L). Data are expressed as 
means ± s.e.m. Asterisks represent p-values as shown: *p < 0.05. 
Statistical significance was measured using a Student's f-test, n = 3 
for each group. 



of the area of the entire apparatus (Dark box insert for mouse 
open field activity, Med Associates Inc). The dark chamber con- 
tains a small arched entrance to allow mice to access the dark 
environment; this entryway was placed toward the center of the 
apparatus so mice could easily access it. Time spent in either the 
light or dark compartment was automatically recorded by Activity 
Monitor software (Med Associates Inc). Relative time spent in the 
light chamber vs. the dark chamber was calculated for each mouse 
by dividing the time spent in either the light or the dark by the 
total exploration time (600 s for all animals). 

ELEVATED PLUS MAZE 

General anxiety behavior was also evaluated using an elevated 
plus maze test as described previously (Chadwick et al, 2011). 
In brief, animals were given 5 min to explore a plus shaped maze 
raised approximately 36 cm off of the ground. Two of the arms 
of the maze were enclosed on all sides; two arms remained open 
for the duration of the test. ANY-Maze Video Tracking Software 
(Stoelting) was used to record time spent in either the closed or 
open arms of the apparatus and the relative time spent in the 
closed or open arms was calculated by dividing the time spent 
in either arm by the total exploration time (300 s for all animals). 

R0TA-R0D TEST 

Motor coordination was evaluated using an accelerating Rota- 
Rod treadmill (Med. Associates, Inc) as described previously 



(Bohlen et al., 2009). Briefly, mice were given two training ses- 
sions on the spinning Rota-Rod apparatus on the day prior to the 
testing day (the first trial was administered in the morning and 
the second trial was administered in the afternoon). Each habitu- 
ation training trial lasted 2 min [4 revolutions per minute (rpm)]. 
On the test day, the mice were placed on the Rota-Rod, which 
gradually accelerated from 4 to 40 rpm over a 5 min time interval. 
The test was performed twice per day and the latency to fall was 
measured and averaged. 

ANIMAL EUTHANIZATION AND TISSUE COLLECTION 

Following completion of behavioral analyses, animals were euth- 
anized via isoflurane overdose (Butler Animal Health Supply, 
Dublin, OH) as described previously and in accord with approved 
animal procedures (Chadwick et al, 201 1). Bodyweight data was 
collected immediately prior to euthanization for each animal. 
Hippocampal and cortical brain tissues were snap frozen on dry 
ice and stored at — 80°C until used for further analyses. 

PROTEIN DIGESTION, iTRAQ LABELING AND WESTERN BLOTTING 

Protein digestion and iTRAQ labeling were processed 
according to the iTRAQ protocol (iTRAQ Reagents — 8plex: 
amine-modifying Labeling Reagents for Multiplexed Relative 
and Absolute Protein Quantitation, ABSciex). All tissue samples 
(Hippocampus and Cortex; WT and BTBR) were prepared in 
parallel throughout the labeling procedures. Briefly, 100 (ig 
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FIGURE 3 | Protein expression in wild type and BTBR mouse 
hippocampus. Western blotting was performed using mouse 
hippocampal protein lysates (A). Protein levels were measured in BTBR 
and control mice with respect to the following proteins in the 
hippocampus: pro-BDNF (B), mature BDNF (C), phospho-TrkB (D), total 



of protein in each condition was acetone precipitated and 
resuspended in 20 |xL of iTRAQ dissolution buffer (0.5 M tri- 
ethylammonium bicarbonate (TEAB), ABSciex) containing 0.1% 
ProteaseMAX detergent (Promega) to denature the proteins. The 
sample was then reduced by adding iTRAQ Reducing Reagent 
[Tris(2-carboxyethyl) phosphine (TCEP), ABSciex] to a final con- 
centration of 5 mM and incubated at 60° C for 1 h. Subsequently, 
the sample was alkylated with iTRAQ Cysteine-Blocking Reagent 
(10 mM methyl methanethiosulfonate (MMTS), ABSciex) for 
lOmin at room temperature. The protein sample was then 
digested with 5 |xg sequencing-grade trypsin (Promega) per 
100 |xg protein at 37°C overnight. Labeling of the samples with 
iTRAQ 8-plex labels was performed at room temperature for 
2 h. After labeling, the samples to be compared were mixed 
and underwent an off-line strong cation exchange (SCX) 
fractionation (ICAT Cation Exchange Buffer Pack, ABSciex) to 
16 fractions to reduce sample complexity. After reversed-phase 
desalting (C18 tips, Pierce), the samples were re-constituted 
in water with 0.1% formic acid, then stored at — 20°C until 
LC/MS/MS analysis. Western blotting procedures as described 
previously (Maudsley et al., 2007; Martin et al., 2009a), were 
performed using the same protein lysates employed for iTRAQ 
labeling. Antibodies used were as follows: rabbit anti-BDNF 
[(brain-derived neurotrophic factor) Santa Cruz Biotechnology, 
Santa Cruz CA] , rabbit antiphosphorylated TrkB and TrkB, rabbit 
anti-phosphorylated Akt and Akt, rabbit antiphosphorylated 



TrkB (E), Akt (F) phospho-synapsin 1 (G), synapsin 1 (H) synaptophysin 
(I), PSD95 (J), spinophilin (K), Neuronal cell marker NeuN (L). Data are 
expressed as means ± s.e.m. Asterisks represent p-values as shown: 
*p < 0.05, **p<0.01. Statistical significance was measured using a 
Student's f-test, n = 3 for each group. 



synapsin 1 and synapsin 1, rabbit anti-spinophilin, rabbit 
anti-PSD95 and rabbit anti-synaptophysin (Cell signaling 
technology, Inc., Danvers, MA), mouse anti-NeuN (EMD 
Millipore Corporation, Billerica, MA) and mouse anti-beta actin 
(Sigma- Aldrich, St. Louis, MO). Results were quantified and 
statistical analysis was performed using an unpaired student 
f-test. p < 0.05 was considered statistically significant. 

LC/MS/MS ANALYSIS 

Samples were analyzed using an Eksigent NanoLC Ultra 2D 
(Dublin, CA) and Thermo Fisher Scientific LTQ Orbitrap XL 
(San Jose, CA). In brief, peptides were first loaded onto a trap 
cartridge (Agilent), then eluted onto a reversed phase PicoFrit 
column (New Objective, Woburn, MA) using a linear 120 min 
gradient of acetonitrile (2-62%) containing 0.1% formic acid 
at 250nL/min flowrate. The eluted peptides were sprayed into 
the LTQ Orbitrap XL. The data-dependent acquisition mode was 
enabled, and each FTMS MSI scan (60,000 resolution) was fol- 
lowed by 6 MS2 scans (alternating CID at unit resolution and 
HCD at 7500 resolution on 3 precursor ions). The spray voltage 
and ion transfer tube temperature were set at 1.8 kV and 180°C, 
respectively. 

DATABASE SEARCH AND iTRAQ QUANTIFICATION 

Proteome Discoverer 1.2 (Thermo Fisher Scientific) was used 
for protein identification and iTRAQ quantification using 
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FIGURE 4 | Quantitative proteomic and bioinformatics analysis of BTBR 
CNS tissues. Log2transformed iTRAQ ratio data (BTBR:WT) snake plot of the 
proteins identified in the cortex (A) and hippocampus (B) are depicted. 
Upregulated proteins are highlighted in red, while downregulated proteins are 



highlighted in green. Three proteins identified by iTRAQ were validated via 
western blot and results correlated to the iTRAQ results (C, D). Data are 
expressed as means ± s.e.m. Asterisks represent p-values as shown: 
*p < 0.05. Statistical significance was measured using a Student's f-test. 



Sequest algorithms. The following criteria were followed: 
SwissProt mouse database; enzyme: trypsin; miscleavages: 2; 
static modifications: methylthio (+45.988 Da on C), iTRAQ8plex 
(+304.205 Da on N-terminus and K); dynamic modifications: 
oxidation (+15.995 Da on M), deamidation (+0.984 Da on N 
and Q); peptide tolerance as 25 ppm; MS2 tolerance as 0.8 Da. 
Peptides reported via all search engines were accepted only if they 
met the false discovery rate of 5%. For iTRAQ quantification, the 
reporter ion intensities in MS2 spectra (m/z 113-121, integration 
width tolerance 50 mmu) were used to calculate the expression 
ratios among the different conditions (Hippocampus and Cortex; 
WT and BTBR). 

BIOINFORMATICS ANALYSIS 

Functional annotational clustering of our proteomic data, i.e., 
gene ontology (GO), KEGG (Kyoto Encyclopedia of Genes and 
Genomes) pathway analysis was performed using WebGestalt 
(http://bioinfo.vanderbilt.edu/webgestalt/, 3/2014) as previously 
described (Chadwick et al., 2010a, 2011). In addition Ingenuity 
Pathway Analysis (http://Ingenuity.com, 3/2014) was employed 
for Canonical Signaling Pathway enrichment investigations. In 
brief, GO analysis allows for broad clustering of genes into 
functional groups, while KEGG pathways analysis allows for 
more specific clustering into signaling pathways (Maudsley 
et al, 2011). Inclusion criteria were set as follows: path- 
way groups needed to meet a minimum population of two 
genes from the input experimental set, and also needed to 



possess a probability significance of enrichment compared to 
a control background dataset of less than p < 0.05 (hyper- 
geometric test of significance). In addition, Venn diagrams 
were also constructed to identify common and uniquely 
altered genes between hippocampus and cortex, using VennPlex 
(Cai et al., 2013), as described previously (Martin et al., 
2009b). 

TEXTROUS! LATENT SEMANTIC ANALYSIS-BASED INVESTIGATION 

Textrous! latent semantic analysis was performed as described 
previously (Chen et al., 2013). In brief, corresponding gene sym- 
bols were uploaded into Textrous! (http://textrous.irp.nia.nih. 
gov/, 3/2014) and then the relevant output indices, i.e., Cosine 
Similarity, Z-score and p-value were extracted. In addition to 
the direct list-based semantic output both the Collective pro- 
cessing mode (generating hierarchical word-clouds), as well the 
Individual processing mode (generating symbol-word heatmaps) 
were applied to provide more nuanced informatic appreciation of 
the specific input datasets. 

STATISTICAL ANALYSIS 

All statistical analyses were conducted using a Student's f-test 
(two-tailed with equal variances: GraphPad Prism, version 5.02). 
p < 0.05 was considered statistically significant throughout the 
study. Error bars represent 95% confidence interval. All data 
represent means ± s.e.m. (standard error of mean). *p < 0.05; 
**p < 0.01; ***p < 0.001. 
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Table 1 | GO term analysis for cortical protein alterations observed in aged BTBR compared to WT control. 
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GO term annotation of the cortical proteins possessing BTBR:WT iTRAQ expression ratios < 1.2 or <0.8 were created using WebGestalt. Each GO term group 
significantly populated (p < 0.05) by at least two independent proteins are depicted. For each GO Term group the specific GO term ID, number of total proteins in 
the curated GO term cluster (C), the number of observed experimental proteins in the GO term group (0), the expected number of proteins in the group based on 
whole-genome background frequency (E), the GO term group enrichment factor {Ft), the GO term enrichment probability (P) and the hybrid score (H * —logm P) sre 
indicated (H). 



RESULTS 

ASD-LIKE BEHAVIORAL PHEN0TYPE PERSISTS IN AGED BTBR MICE 

Commensurate with an extant ASD phenotype aged BTBR mice 
demonstrated no sociability preference between an object or 
another mouse, while WT controls spent significantly more time 
with another mouse compared to an object (Figure 1A). Relative 
percentage of time spent with either the mouse or the object 
was used to calculate a social discrimination index [DI: (novel 
time/total time-same time/total time) x 100] (Figure 2B). WT 
and BTBR mice performed similarly in the novel object prefer- 
ence test of learning and memory, with both groups spending 
significantly more time with a novel object compared to a familiar 



object (Figures 1C,D). Both WT and BTBR mice significantly 
preferred to explore the peripheral areas of the chamber over 
the center of the chamber, though the BTBR mice did show a 
significantly higher preference for the periphery over the center 
(Figure IE). No significant differences in total distance trav- 
eled (Figure IF) or number of ambulatory episodes (Figure 1G) 
were found. BTBR mice however demonstrated significantly less 
jumping activity compared to WT (Figure 1H). In light-dark 
exploration tests BTBR mice spent significantly more time in the 
dark chamber compared to the lighted chamber, while WT con- 
trol mice show no statistical difference in preference for either 
chamber (Figure II). In the elevated plus maze test of anxiety, 
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Table 2 | GO term analysis for hippocampal protein alterations observed in aged BTBR compared to WT controls. 
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GO Class 


GO Term 


GO ID 


C 


E 


O 


R 


P 


H 


Cellular component 


Neuron projection 




DO 1 


JO 


7 OR 
/ .ZD 


A RR 
H. DO 


1 nop 11 
I .uzt- I I 


rci ninQ7 

OU.U I Uo / 


Cellular component 


Synaptic vesicle 


QCiO 1 
OUZ I 


1 nc; 
lUo 


1 n 
I u 


1 17 
I.I/ 


Q RR 
O.DD 


o 1 a p 

Z. I'+E— uo 


AQ /17/IQK 

4o.4 /^-yo 


Cellular component 


Axon 


on a o a 

30424 


286 


19 


3.19 


5.96 


o.4/b-uy 


4o.1Uyo1 


cellular component 


Respiratory chain 


70469 


76 


8 


0.85 


9.45 


-i oo n r\ o 

1 .Jzb-Oo 


nn unco 

4o.nUoo 


Cellular component 


Clathrin-coated vesicle membrane 


onccc; 
oUbbD 


1 1 n 
I I U 


1 n 
IU 


1 OQ 

I .Zo 


o. I O 


q o 1 p nc 
o.Z 1 1— U0 


A A QOROO 

44.oZ0yZ 


Cellular component 


Vesicle membrane 


1 2506 


355 


2 1 


3.96 


5.31 


o a 7 p no 


a o ocon A 
4Z.ooZy4 


Cellular component 


Cytoplasmic vesicle membrane 


oUboy 


OAO 

o4Z 


on 
ZU 


Q Q1 
O.O I 


R OR 
O.ZO 


o coc no 
Z.DZt— Uo 


on qoocc 

jy .oyzoo 


Cellular component 


Cytoplasm 


5737 


91 30 


1 70 


101 .72 


1 .67 


1 .52E-24 


OO 77COO 

oy. / /ooZ 


Cellular component 


Cytoplasmic part 


A A A A A 

44444 


6772 


143 


75.45 


1 .9 


1 OC P 0 1 

I .oob-Z I 


On 0Q700 


Cellular component 


Clathrin-coated vesicle 


oU I So 


207 


14 


2.31 


6.07 


7.60E-07 


07 1 yl O A R 

o /. 1 4o4b 


Cellular component 


Cell projection 


a ooqc; 

4zyyo 


1 oon 
I ZoU 


A R 

40 


1*3 "7 
lo. / 


Q OO 

o.Zo 


O R 1 P 11 
Z.D I t— I I 


OA 7CQH7 

o4. /oyu/ 


Cellular component 


Coated vesicle membrane 


oncco 
JUbbz 


1 50 


1 1 


1 .67 


6.58 


c 70.P nc 
o. /yb-Uo 


o/i nnco 
34.UUb3 


Cellular component 


Cytoplasmic vesicle part 


AAA OO 

444oo 


401 


2 1 


4.47 


4.7 


c ice no 
O. I Ob— Uo 


OO 0QQ07 


Cellular component 


Cytosol 


5829 


2372 


69 


26.43 


2.61 


4.11 E-13 


32.32787 


Cellular component 


Coated vesicle 


30135 


254 


15 


2.83 


5.3 


1 .44E-06 


30.96068 


Cellular component 


Mitochondrion 


5739 


1525 


48 


16.99 


2.83 


5.17E-10 


26.28082 


Cellular component 


Synapse part 


44456 


370 


18 


4.12 


4.37 


1 .46E-06 


25.50178 


Cellular component 


Membrane-bounded vesicle 


31988 


883 


32 


9.84 


3.25 


4.82 E-08 


23.7801 


Cellular component 


Mitochondrial inner membrane 


5743 


354 


17 


3.94 


4.31 


3.58E-06 


23.47276 


Cellular component 


Cytoplasmic membrane-bounded 
vesicle 


16023 


862 


31 


9.6 


3.23 


8. 13 E-08 


22.90041 


Cellular component 


Mitochondrial part 


A A aoo 
444zy 


"7/1 A 

/44 


o~7 
Z / 


Q OG 

o.zy 


o.Zo 


R OOP n7 
D.oZt— U/ 


on oocor; 
ZU.oZboo 


Cellular component 


Organelle inner membrane 


1 flOCC 


384 


1 7 


4.28 


3.97 


1 m p nc 
I .U I L— Uo 


1 n oooo a 
iy.oJZo4 


Cellular component 


Vesicle 


Q 1 OQO 

o i yoz 


y / u 


QO 
oZ 


1 n qi 
I U.o I 


O OR 

z.yo 


q onp m 
o.ZUt— U / 


1 Q O O /I 7C 

i y.zz4/b 


Cellular component 


Cytoplasmic vesicle 


31410 


928 


31 


10.34 


3 


o nop n7 
o.yoL-U/ 


1 n o 1 coo 
I y.Z I boZ 


Cellular component 


Intracellular part 


A A Al A 

444Z4 


1 2237 


1 79 


136.34 


1 .31 


O A O P 1 O 

Z.4ZL- lo 


1 C C070 

Ib.OZ /Z 


Cellular component 


Intracellular 


5622 


1 2564 


1 80 


i on no 


1 .29 


1 OOP 1 o 

I .oZL- Iz 


1 c OO A A C 

1 0.oZ44b 


Cellular component 


Organelle membrane 


o 1 nnn 

o iuyu 


2373 


57 


26.44 


2.16 


o 1 o p no 
o. I JL- Uo 


1 rr o-| yio 
! O.O I4Z 


Cellular component 


Intracellular organelle part 


44446 


6725 


120 


74.93 


1 .6 


b.oUb-IU 


14. / / ooz 


Cellular component 


Organelle part 


444 zz 


6812 


120 


75.9 


1 .58 


1 .4zb-Uy 


13.9 /yJo 


Cellular component 


Cytoskeleton 


5856 


1790 


45 


19.94 


2.26 


1 .Uob-Uo 


1 O (7000-1 

lo.oUZol 


Cellular component 


Protein complex 


43234 


3278 


70 


36.52 


1.92 


9.89E-08 


13.44922 


Cellular component 


Cell part 


A A A r* A 

44464 


14643 


189 


163.15 


1 .16 


-i nAr -in 

1 .90L-10 


■1 -1 O — n f— 

11.27665 


Cellular component 


Cell 


5623 


14644 


189 


163.16 


1.16 


1.90E-10 


11.27665 


Cellular component 


Macromolecular complex 


32991 


3864 


76 


43.05 


1.77 


5.44E-07 


11.08799 


Cellular component 


Intracellular organelle 


43229 


10636 


155 


118.5 


1.31 


8.09E-08 


9.290587 


Cellular component 


Organelle 


43226 


10651 


155 


118.67 


1.31 


8.13E-08 


9.287781 


Cellular component 


Intracellular membrane-bounded 
organelle 


43231 


9587 


138 


106.82 


1.29 


1.32E-05 


6.29446 


Cellular component 


Membrane-bounded organelle 


43227 


9598 


138 


106.94 


1.29 


1 .36E-05 


6.277735 



GO term annotation of the cortical proteins possessing BTBR:WT iTRAQ expression ratios < 1.2 or <0.8 were created using WebGestalt. Each GO term group 
significantly populated (p<0.05) by at least two independent proteins are depicted. For each GO Term group the specific GO term ID, number of total proteins in 
the curated GO term cluster (C), the number of observed experimental proteins in the GO term group (0), the expected number of proteins in the group based on 
whole-genome background frequency (E), the GO term group enrichment factor (Fi), the GO term enrichment probability (P) and the hybrid score (R * —logw P) 3re 
indicated (H). 



both WT and BTBR mice significantly preferred closed arms 
to open arms (Figure 1J). Motor coordination, assessed using a 
Rotarod test was similar in both groups (Figure IK). 

ALTERED PROTEIN EXPRESSION IN CNS OF AGED BTBR MICE 

Cortical and hippocampal western blotting of proteins, selected 
for their ability to indicate alterations of neurosynaptic activity, 



was employed to detect potential expression differences between 
aged WT and BTBR mice. Our panel of neurosynaptic factors 
analyzed in cortical tissue is indicated in Figure 2A. No signifi- 
cant difference between WT and BTBR pro-BDNF (brain-derived 
neurotrophic factor) levels was detected (Figure 2B), whereas a 
significant decreased in BTBR mature BDNF was observed com- 
pared to WT (Figure 2C). Levels of phosphorylated TrkB receptor 
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FIGURE 5 1 Signaling pathway analysis of cortical and 
hippocampal proteins differentially regulated BTBR and WT. 

Cortical and hippocampal proteins exhibiting differential regulation 
between BTBR and control mice were analyzed for potential 
functional pathway interactions using KEGG (A-cortex, 



B-hippocampus) and Ingenuity Canonical Signaling pathway analysis 
(C-cortex, D-hippocampus). For (A,B), the top 10 most 
significantly-populated (calculated with Hybrid scores: enrichment 
factor (R) * (— logio enrichment probability, P) neuronally-specific 
KEGG or Canonical signaling pathways are shown. 



were similar in WT and BTBR cortex (Figure 2D) while the total 
TrkB expression was significantly lower in BTBR mice compared 
to WT (Figure 2E). No significant differences in Aktl expres- 
sion was noted in the BTBR model compared to WT (Figure 2F). 
Levels of phosphorylated synapsin 1 (Figure 2G), but not total 
synapsin 1 (Figure 2H) were significantly elevated in BTBR mice 
compared to WT. Expression of both synaptophysin (Figure 21) 
and PSD95 (Figure 2 J) was not affected in BTBR mice, while 
both spinophilin (Figure 2K) and NeuN (Figure 2L) levels were 
significantly elevated in BTBR mice compared to WT. 

In the hippocampus we again found no significant differences 
in the levels of pro-BDNF between BTBR and WT (Figure 3B) 
that was coincident with a significant reduction in extant 
BDNF levels (Figure 3C). Phosphorylated TrkB expression was 
unchanged in BTBR mice compared to WT (Figure 3D), while 
there was a trend, similar to that in the cortex, to reductions of 
total TrkB receptor expression in the BTBR mice (Figure 3E). As 
with the cortex no significant alterations in Aktl were observed 
(Figure 3F). Significant potentiation of both phosphorylated 
(Figure 3G) and non-phosphorylated synapsin 1 (Figure 3H) 
expression levels was seen in the BTBR hippocampus. In 
accordance with the cortical data we found no significant 
BTBR-induced changes in the hippocampal expression of synap- 
tophysin (Figure 31) or PSD95 (Figure 3J). Levels of both 



spinophilin (Figure 3K) and NeuN (Figure 3L) were again sig- 
nificantly potentiated in the BTBR mice, mirroring our previous 
cortical data. 

SYSTEMIC ANALYSIS OF BTBR CENTRAL NERVOUS PROTEIN 
EXPRESSION 

To gain an unbiased appreciation of the multiple tissue pro- 
tein expression pattern changes generated in the BTBR mice 
compare to WT we applied quantitative isobaric mass-tag label- 
ing (iTRAQ) to cortical and hippocampal tissue extracts. In the 
cortex we identified and generated relative expression profiles 
for 674 proteins (Figure 4A: Table SI), while in the hippocam- 
pus 656 proteins were identified and quantified in BTBR mice 
(Figure 4B: Table S2). Log2 -transformed iTRAQ expression ratios 
were employed to generate snake plot graphs of individual pro- 
teins. To initially validate some of the iTRAQ data we chose three 
random proteins, Rab3A, CaMKIID and Cplxl and found that at 
the western level of detection our results were comparable to our 
iTRAQ expression ratio data (Figures 4C,D). 

GENE ONTOLOGY BIOINFORMATIC ANALYSIS OF BTBR-SPECIFIC 
ALTERED PROTEINS 

Using a standard arbitrary cut-off of >1.2 (elevated in BTBR 
vs. WT) or <0.8 (reduced in BTBR vs. WT) for the iTRAQ 
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Table 3 | KEGG Pathway analysis for cortical proteins altered in BTBR mice compared to WT controls. 



KEGG pathway 


C 


E 


O 


Ft 


P 


H 


Glutathione metabolism 


50 


2 


0.07 


27.38 


0.0188 


47.25356 


MTOR signaling pathway 


52 


2 


0.08 


26.33 


0.0188 


45.44142 


Arginine and proline metabolism 


54 


2 


0.08 


25.35 


0.0188 


43.7501 


Shigellosis 


61 


2 


0.09 


22.44 


0.0188 


38.7279 


Renal cell carcinoma 


70 


2 


0.1 


19.56 


0.0188 


33.75747 


Chronic myeloid leukemia 


73 


2 


0.11 


18.75 


0.0188 


32.35954 


Progesterone-mediated oocyte maturation 


86 


2 


0.13 


15.92 


0.0188 


27.47541 


Gap junction 


90 


2 


0.13 


15.21 


0.0188 


26.25006 


Melanogenesis 


101 


2 


0.15 


13.56 


0.0192 


23.27844 


GnRH signaling pathway 


101 


2 


0.15 


13.56 


0.0192 


23.27844 


Oocyte meiosis 


112 


2 


0.16 


12.22 


0.0192 


20.97806 


Leukocyte transendothelial migration 


116 


2 


0.17 


11.8 


0.0192 


20.25705 


Vascular smooth muscle contraction 


116 


2 


0.17 


11.8 


0.0192 


20.25705 


Neurotrophin signaling pathway 


127 


2 


0.19 


10.78 


0.0192 


18.50601 


Spliceosome 


127 


2 


0.19 


10.78 


0.0192 


18.50601 


Tight junction 


132 


2 


0.19 


10.37 


0.0192 


17.80217 


Hepatitis C 


134 


2 


0.2 


10.22 


0.0192 


17.54466 


Natural killer cell mediated cytotoxicity 


136 


2 


0.2 


10.07 


0.0192 


17.28716 


Insulin signaling pathway 


138 


2 


0.2 


9.92 


0.0192 


17.02965 


Alzheimer's disease 


167 


2 


0.24 


8.2 


0.0261 


12.98355 


Chemokine signaling pathway 


189 


2 


0.28 


7.24 


0.0313 


10.89226 


Metabolic pathways 


1130 


6 


1.65 


3.63 


0.0188 


6.264807 



KEGG signaling pathway annotation of the cortical proteins possessing BTBR:WT iTRAQ expression ratios < 1.2 or < 0.8 was created using WebGestalt. Each KEGG 
pathway significantly populated (p < 0. 05) by at least two independent proteins are depicted. For each populated KEGG pathway the number of total proteins in the 
curated KEGG pathway (C), the number of observed experimental proteins in the KEGG pathway 10), the expected number of proteins in the pathway based on 
whole-genome background frequency (E), the KEGG pathway enrichment factor (R), the KEGG pathway enrichment probability (P) and the hybrid score (R * —logw 
P) are indicated (H). 



ratios for differentially expressed proteins between BTBR and 
WT mice we investigated the functional clustering of the BTBR- 
specific protein changes into Gene Ontology term groups. 
Analysis of the differentially-regulated [elevated (<1.2) and 
reduced (<0.8)] proteins in the BTBR cortex we found a 
strong GO-biological process clustering of endoplasmic reticu- 
lar, vesicular transport and neurotransmission related functions 
(Table 1). No GO-molecular function annotation was achieved 
with the cortical data however. With respect to GO-cellular 
compartment analysis we found the differential BTBR proteins 
were strongly associated with endocytic activity, neuronal out- 
growth and architecture as well as mitochondrial metabolic 
activity (Table 1). GO-biological process analysis of the dif- 
ferential BTBR hippocampal protein set indicated a strong 
bias toward neuronal synaptic secretion, neuronal plasticity, 
metabolic activity and inter-cellular contact (Table 2). GO- 
cellular compartment interrogation of the hippocampal dataset 
revealed a strong bias for protein interactions in the synap- 
tic cleft, clathrin-coated vesicles and mitochondrial-associated 
metabolic areas (Table 2). GO-molecular function annotation 
was evident with the hippocampal data and revealed the 
presence of GO groups associated with malate metabolism, 
redox activity, synaptic vesicle docking and structural modeling 
(Table 2). 



SIGNALING PATHWAY ANALYSIS OF DIFFERENTIALLY REGULATED 
BTBR-SPECIFIC PROTEINS 

While GO annotation indicates some degree of functional com- 
monality between groups of proteins, associating multiple pro- 
teins with classically-identified molecular signaling paradigms 
provides additional functional information regarding the poten- 
tial physiological ramifications of specific group of differentially- 
regulated factors. We therefore performed both KEGG and 
Ingenuity-base pathway analysis upon our cortical and hip- 
pocampal BTBR protein sets. Inspecting the top ten highest scor- 
ing KEGG pathways from cortical proteins (Figure 5 A: Table 3) 
a strong population of metabolic (Glutathione metabolism, 
Insulin signaling pathway, Metabolic pathways), neurotransmis- 
sion/degenerative (Neurotrophic signaling pathway, Alzheimer's 
disease) and ultrastructural (Gap junction, Tight junction, 
Spliceosome) signaling pathways was evident. Applying the same 
analysis to the hippocampus an intense population of signaling 
systems linked to oxidative metabolism [Oxidative phosphoryla- 
tion, Pyruvate metabolism, Citrate cycle (TCA cycle)] and cen- 
tral neurodegeneration (Alzheimer's disease, Huntington's disease, 
Parkinson's disease, Long-term potentiation) was clear (Figure 5B: 
Table 4). Extraction of potential physiological meaning from 
complex data sets is best performed using multiple informatics 
tools. Therefore, we also sought interpretation of the potential 
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Table 4 | KEGG Pathway analysis for hippocampal proteins altered in BTBR mice compared to WT controls. 



KEGG pathway 


C 


E 


O 


R 


P 


H 


Oxidative phosphorylation 


132 


1 5 


0.61 


24.63 


1 .42 E— 1 5 


occ con 1 

job.oyy i 


Alzheimer's disease 


167 


1 7 


0.77 


22.06 


1 no c 1c 
I .UJt— I O 


QCQ C"7CQ 

jDz.b/bo 


Huntington's disease 


1 QQ 
I O J 


1 Q 
I O 


n Q/i 
U.o4 


O 1 QQ 

Z I . JZ 


Q QQ C 1 "7 

z.yz t— i / 


JdZ.D I O 


Parkinson's disease 


1 on 
I JU 


1 Q 

l Z 


n c 
U.b 


on 
ZU 


1 QQ P 11 
I .oot— I I 


Q 1 /I CO 

Z I4.D loo 


Long— term potentiation 


70 


8 


0.32 


24.77 


1 "7/i C" no 
1 . /4t-Uo 


-1 no om c 

i yz.zu I b 


Gastric acid secretion 


~IA 


Q 
O 


n Q/i 
U. J4 


QQ A Q 
ZJ.4J 


q Q/i p no 
Z. J4t— Uo 


1"7Q "7QQQ 

i /o. /oyj 


Collecting duct acid secretion 


07 
Z I 


A 

4 


n n 
U. IZ 


QQ 1 1 
JZ. I I 


4.U I t— Uo 


1/111 QQ 
14 I. I oJ 


Pyruvate metabolism 


40 


5 


0.18 


27.09 


7.57E-06 


-IOO ~70 (T O 

1 Jo. /zdJ 


Metabolic pathways 


1130 


37 


5.21 


7.1 


C QQ C" 10 

0. Jot— i y 


1QO "711/1 

1 zy. / 1 14 


GnRH signaling pathway 


101 


8 


0.47 


17.17 


z.1yt-U/ 


11/1 O /I /I c 

114.J44b 


Oocyte meiosis 


1 1 2 


8 


0.52 


1 5.48 


a r\A p n~7 
4.U4t— U / 


QQ 0"7QQO 

yy.y / jzz 


Calcium signaling pathway 


177 


10 


0.82 


1 2.24 


1 1 /i c n~7 
I . I4£— U/ 


OA QQO A O 

o4.yoJ4o 


Vasopressin-regulated water reabsorption 


44 


4 


0.2 


19.7 


U.UUUz 


~70 opmi 

/z.oby / I 


Endocytosis 


201 


10 


0.93 


10.78 


q no c n~7 
J.UJt— U / 


"7n oinnc 
/U.Z /UUd 


Citrate cycle (TCA cycle) 


30 


3 


0.14 


21 .67 


0.001 


65.01 


Butanoate metabolism 


on 
JU 


J 


n 1 a 
U. 14 


Z 1 .0/ 


n nm 
U.UU I 


bD.U I 


Propanoate metabolism 


32 


3 


0.15 


20.32 


0.001 


60.96 


Alanine, aspartate and glutamate metabolism 


JZ 


Q 
O 


U. lb 


on QQ 
ZU. JZ 


n nm 
U.UU I 


cn nc 

bu.yb 


Melanogenesis 


IU1 


6 


0.47 


1 2.87 


4.1 5E— 05 


QQC"7Q 

ob. jyo/z 


Vibrio cholerae infection 


54 


4 


0.25 


1 6.05 


n nnn/i 
U.UUU4 


n /i coco/1 
D4.bJby4 


Shigellosis 


p. 1 


4 


n qq 

U.Zo 


1 /l Q 1 
I4.Z I 


n nnn7 
u.uuu / 


A A QQ1 1 K 
44.0 J I I b 


Wnt signaling pathway 


1 oU 


-7 

I 


n rq 

u.oy 


1 n 1 1 
IU. I I 


4.U 1 1— Uo 


/I /I /I KOO 1 

44.4bZZ I 


Sa ivary secretion 


89 


5 


0.41 


12.18 


n nnno 
U.UUUJ 


a o nnocc 
4z.yUobb 


Gap junction 


on 


O 


n /i 0 
U.4Z 


1 q n/i 
1 Z.U4 


n nnnQ 
U.UUUJ 


AO A"\ K/l C 

4Z.4 I b4b 


Beta-Aanine metabolism 


22 


2 


0.1 


1 9.7 


0.0071 


/i o oono 1 
4z. JJUz I 


Glioma 


65 


4 


0.3 


13.34 


0.0007 


/i 0 nn^on 

4z.UobJy 


Neurotrophin signaling pathway 


1 27 


6 


0.59 


10.24 


U.UUU1 


40.96 


Valine, leucine and isoleucine degradation 


44 


3 


0.2 


14.78 


U.UUzJ 


00 nno/^r 1 

Jo.9yJbb 


Proteasome 


44 


3 


0.2 


14.78 


n nnQQ 
U.UUzJ 


OO OQOCC 

jo.yyjbb 


Epithelial cell signaling in Helicobacter pylori infection 


68 


4 


0.31 


12.75 


n nnno 

u.uuuy 


OO OOQ /I 1 
Jo.OJJ4 I 


PPAR signaling pathway 


70 


4 


0.32 


1 2.38 


n nnno 

u.uuuy 


Q"7 "7nC/1 Q 

J /. /Ub4o 


Long-term depression 


/u 


A 
4 


n qq 

U. JZ 


1 9 QQ 
I Z.Jo 


n nnna 
u.uuuy 


Q7 7nR/l Q 
J /. / Ub4o 


Bile secretion 


"71 
/ I 


A 

4 


n qq 
U.Jo 


1 Q Q 1 
I Z.Z I 


n nnnQ 

u.uuuy 


Q7 -1 QQ"7 

J /. 1 00 / 


Cardiac muscle contraction 


7-7 


A 
4 


n oc 
U. Jo 


I I .Zo 


n nm o 
U.UU I Z 


QQ QQQ/1 Q 
JZ.0004Z 


Phosphatidylinositol signaling system 


7R 
/o 


/I 
4 


n qr 
U. Jo 


11 11 
I I . I I 


n nm o 

U.UU I z 


QQ A P^QQ 
JZ .4DUJ 


Pentose phosphate pathway 


27 


2 


0.12 


1 6.05 


n m no 
U.U IUJ 


OI OQOQC 

j i.oyjyb 


Amyotrophic lateral sclerosis (ALS) 


53 


3 


0.24 


1 2.27 


n nnoc 
U.UUJb 


QO no A 1 7 

zy.yo4 1 / 


ErbB signaling pathway 


87 


4 


0.4 


9.96 


n nm c 
U.UUlo 


Q7 Q/i cnc 
z /.o4byb 


Inositol phosphate metabolism 


O / 


Q 
J 


n or 

U.ZO 


11/1 
1 1 .4 1 


n nn/iQ 
U.UU4J 


Q7nno 1 q 

Z /.UUZ I z 


Regulation of actin cytoskeleton 


213 


7 


0.98 


7. 1 2 


n nnno 
U.UUUJ 


on noon 
zb.UoZy 


Glycolysis/Gluconeogenesis 


65 


3 


0.3 


10 


0.0057 


oo /i/i nc 
ZZ. 44 Izb 


Prion diseases 


Of: 


0 
Z 


u. lb 


1 Q QQ 
I Z.Jo 


U.U I DO 


qq onncQ 
ZZ. JUUbJ 


African trypanosomiasis 


oc 
oD 


0 
Z 


n 1 c 
U. I 0 


1 Q QQ 
I Z.Jo 


U.U I DO 


qq QnncQ 
ZZ. JUUbJ 


Pancreatic secretion 


1 n-i 
IU1 


A 

4 


0.47 


8.58 


n nnoc 
U.UUzo 


QQ 17QCO 

zz. I /ybJ 


Adipocytokine signaling pathway 


68 


3 


0.31 


9.56 


n nnc/i 
U.UUo4 


on 07000 

zu.y /zyz 


Bacterial invasion of epithelial cells 


70 


3 


0.32 


9.29 


0.0066 


20.25644 


Renal cell carcinoma 


70 


3 


0.32 


9.29 


0.0066 


20.25644 


Vascular smooth muscle contraction 


116 


4 


0.54 


7.47 


0.0039 


1799475 


Antigen processing and presentation 


76 


3 


0.35 


8.55 


0.0081 


17.88245 


Purine metabolism 


162 


5 


0.75 


6.69 


0.0023 


17.65004 


Protein processing in endoplasmic reticulum 


165 


5 


0.76 


6.57 


0.0023 


17.33345 


Lysosome 


121 


4 


0.56 


7.16 


0.0044 


16.87288 
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KEGG pathway 
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E 


O 


R 


P 


H 


Spliceosome 


1 27 


4 


0.59 


6.83 


0.005 


1 c "7 1 cno 
I b. / I bUo 


Chemokine signaling pathway 


189 


5 


0.87 


5.73 


U.UUob 


I4.UUZ09 


Rheumatoid arthritis 


y i 




n a o 
U.4Z 


~7 1 A 


U.U I Z / 


Io.booo4 


Fc gamma R-mediated phagocytosis 


94 


3 


0.43 


6.92 


U.U loo 


n 01 cn 1 

i z.y i by i 


TastG transduction 


pro 
DZ 


Z 


n oa 
U.Z4 


Q OA 

o.o4 


U.Uo I o 


1 O K/1 "71 ft 

I z.b4/ I b 


MAPK signaling pathway 


Zoo 


ft 
0 


1 OA 

I.Z4 


A 

4.0D 


U.UUoz 


^o inrim 
I Z. IUUUZ 


Arginine and proline metabolism 


54 


2 


0.25 


8.03 


0.0325 


11.94958 


Pathogenic Escherichia coli infection 


56 


2 


0.26 


7.74 


0.0342 


11.34664 


Chagas disease (American trypanosomiasis) 


104 


3 


0.48 


6.25 


0.017 


11.05969 


Amoebiasis 


106 


3 


0.49 


6.13 


0.0177 


10.73992 


Toxoplasmosis 


132 


3 


0.61 


4.93 


0.0309 


7444505 


Natural killer cell mediated cytotoxicity 


136 


3 


0.63 


4.78 


0.0324 


7.119595 


Phagosome 


153 


3 


0.71 


4.25 


0.0416 


5.868853 



KEGG signaling pathway annotation of the hippocampai proteins possessing BTBR:WT iTRAQ expression ratios < 1.2 or < 0.8 was created using WebGestalt. Each 
KEGG pathway significantly populated (p < 0.05) by at least two independent proteins are depicted. For each populated KEGG pathway the number of total proteins 
in the curated KEGG pathway (C), the number of observed experimental proteins in the KEGG pathway (O), the expected number of proteins in the pathway based 
on whole-genome background frequency (E), the KEGG pathway enrichment factor (R), the KEGG pathway enrichment probability (P) and the hybrid score (R * 
-logw P) are indicated (H). 



signaling pathways enriched in the BTBR tissues using canon- 
ical signaling pathway analysis (Ingenuity). Canonical signaling 
analysis of the cortex revealed a strong bias toward cytokine 
signaling (Leptin signaling in obesity, JAK family kinases in 
IL-6-type signaling), neurodegenerative disease related activity 
(Amyloid processing, CDK5 signaling, CNTF signaling) and energy 
metabolism (IGF-1 signaling) (Figure 5C: Table 5). Similar path- 
way processing of hippocampai data again revealed a strong 
prediction of metabolism-related activity (Oxidative phospho- 
rylation, Mitochondrial dysfunction, G protein signaling medi- 
ated by Tubby), neurodegenerative activity (Huntington's disease 
signaling) and interestingly, given the strong presentation of, 
Melanocyte development and pigmentation signaling' in the cor- 
tex (Figure 5C), melatonin-related activity (Melatonin signaling) 
(Figure 5D: Table 6). 

BTBR-REGULATED COHERENT PROTEIN REGULATION PATTERNS 

Along with our global analysis of protein set data from the two 
tissues we also investigated the predicted functional nature of the 
divergent expression polarity subsets of BTBR-specific proteins 
uniquely regulated in either the cortex or hippocampus. Using 
VennPlex (Figure 6A: Table 7) (Cai et al., 2013) we were able to 
identify multiple groups of coherently-controlled proteins across 
multiple CNS tissues. Using our previously developed informat- 
ics application, Textrousl, we generated physiological predictions 
extracted from multiple databases using latent semantic analysis 
(Chen et al., 2013). Using the collective processing capacity of 
Textrousl, which generates interactive hierarchical word-clouds, 
we found that the tissue-unique up or down regulated sets of pro- 
teins generated quite distinct functional signatures in the cortex 
or hippocampus. For example, the specifically upregulated pro- 
teins in the cortex were strongly focused into activity controlling 
excitatory amino acid synaptic activity (Figure S1A: Table S3), 
while the downregulated cortical only proteins demonstrated a 



strong link to a reduction of protein kinase activity (Figure SIB: 
Table S4). With respect to the hippocampus, the upregulated pro- 
tein set appeared to be linked to increases in heat-shock factors 
and synaptic pathophysiology (Figure S2A: Table S5), while the 
downregulated hippocampai protein set suggested a reduction in 
oligodendrocyte generation and CNS myelination (Figure S2B: 
Table S6). 

In addition to our global data set informatics interpreta- 
tion we also sought to investigate the function of coherently- 
regulated groups of BTBR-associated proteins observed across 
multiple tissues. We therefore focused on the subset of proteins 
(12 upregulated; 4 downregulated: Table 7) coherently altered 
in a BTBR-specific manner in both the cortex and hippocam- 
pus (Figure 6B). To explore the potential functional connec- 
tion between these select proteins we applied our novel infor- 
matic platform Textrousl (Chen et al, 2013). Textrousl allows 
the extraction of functional data, using latent semantic anal- 
ysis (LSA), from even the smallest of datasets. Textrousl is 
a unique LSA tool as it possesses two forms of functional 
data output, i.e., collective and individual processing. Collective 
processing of the coherent BTBR protein set results in the 
creation of a hierarchical word-cloud revealing a strong and 
novel link between the BTBR ASD phenotype and N-WASP 
(neuronal Wiskott-Aldrich Syndrome protein) mediated activ- 
ity (Figure 6C: Table 8). In addition there was a second, less 
strong demonstration of the importance of inherited genomic 
activity in ASD phenotypes. The individual processing mode 
in Textrousl also investigates the strongest individual protein- 
term relationships within a specific protein subset (Figure 6D). 
Again with this distinct analytical output we observed a strong 
relationship between cytoskeletal remodeling activity (f-actin, 
scaffolding, n-wasp, filopodia) and proteins such as Itsnl (inter- 
sectin 1), Ap3b2 (adaptor-related protein complex 3, beta 
2 subunit), Dbnl (drebrinlike) and Ermn (Ermin, ERM-like 
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Table 5 | Ingenuity Canonical Signaling Pathway analysis of proteins 
differentially regulated in the cortex of aged BTBR compared to WT 
controls. 

Canonical pathways log(p-value) Ratio 10 x Hybrid 
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Canonical pathways 


loglp- value) 


Ratio 


10 x Hybrid 


Relaxin Signaling 


1.98E - 


- 00 


1.83E- 


-02 


0.36234 


FGF Signaling 


1.49E - 


- 00 


2.13E- 


-02 


0.31737 


Gap Junction Signaling 


1.83E - 


- 00 


1 66 E- 


-02 


0.30378 


Apoptosis Signaling 


1.47E - 


- 00 


9 nnF 


n? 
-uz 


0.294 


Acute Phase Response 


1.74E H 


- 00 


I.DDC- 


09 

-uz 


0.28884 


Signaling 












G-Protein Coupled Receptor 


1.99E - 


- 00 


1 45E- 


-02 


0.28855 


Signaling 












Mouse Embryonic Stem Cell 


1.41 E - 


- 00 


2.02E- 


-02 


0.28482 


Pluripotency 












Neuregulin Signaling 


1.47E + 00 


1.92E- 


-02 


0.28224 


FAK Signaling 


1.48E - 


- 00 


1.89E- 


-02 


0.27972 


Chronic Myeloid Leukemia 


1.42E - 


- 00 


1.89E- 


-02 


0.26838 


Signaling 












VEGF Signaling 


1.43E - 


- 00 


1.83E- 


-02 


0.26169 


IL-1 Signaling 


1.42E - 


- 00 


1.83E- 


-02 


0.25986 


Sertoli Cell-Sertoli Cell 


1.69E - 


- 00 


1.52E- 


-02 


0.25688 


Junction Signaling 












PPARa/RXRa Activation 


1.69E - 


- 00 


1.50E- 


-02 


0.2535 


HGF Signaling 


1.37E - 


- 00 


1.80E- 


-02 


0.2466 


Mitochondrial Dysfunction 


1.76E H 


- 00 


1.40E- 


-02 


0.2464 


Ephrin Receptor Signaling 


1.70E - 


- 00 


1.43E- 


-02 


0.2431 


Role of NFAT in Cardiac 


1.66E - 


- 00 


1.44E- 


-02 


0.23904 


Hypertrophy 












Calcium Signaling 


1.67E - 


- 00 


1.38E- 


-02 


0.23046 


Integrin Signaling 


1.57E - 


- 00 


1.44E-02 


0.22608 


Leukocyte Extravasation 


1.56E - 


- 00 


1.43E-02 


0.22308 


Signaling 












Breast Cancer Regulation by 


1.59E - 


- 00 


1.40E- 


-02 


0.2226 


Stathminl 












Fc Epsilon Rl Signaling 


1.30E - 


- 00 


1.71 E- 


-02 


0.2223 


Nitric Oxide Signaling in the 


1.36E - 


- 00 


1.60E- 


-02 


0.2176 


Cardiovascular System 












NGF Signaling 


1.32E - 


- 00 


1.64E- 


-02 


0.21648 


Axonal Guidance Signaling 


1.86E - 


- 00 


1.03E- 


-02 


0.19158 


Cardiac Hypertrophy 


1.42E - 


- 00 


1.20E- 


-02 


U. ! /U4 


Signaling 












Molecular Mechanisms of 


1.58E - 


- 00 


1.03E- 


-02 


0.16274 


Cancer 












Colorectal Cancer Metastasis 


1.34E - 


- 00 


1.12E- 


-02 


0.15008 


Signaling 












Protein Kinase A Signaling 


1.43E - 


- 00 


9.78E- 


-03 


0.139854 



Canonical signaling pathway annotation of the cortical proteins possessing 
BTBR:WT iTRAQ expression ratios <12 or <0.8 was created using Ingenuity 
Pathway Analysis. Each signaling pathway significantly populated (p < 0.05) by 
at least two independent proteins are depicted. For each populated signaling 
pathway the—logw of the enrichment probability (p-value), the enrichment ratio 
and the calculated hybrid score [10 * (enrichment ratio * —logw(p-value))] are 
depicted. 

protein) with the BTBR phenotype. Interestingly other tight 
relationships between BTBR coherently-regulated factors were 
found, e.g., myelin-associated activity [Lancll (LanC bacterial 
lantibiotic synthetase component C)-like 1 and Ermn] and 
inherited neuronal atrophy and impairment [Prrt2 (proline-rich 



Leptin Signaling in Obesity 


3.97E - 


F 00 


4.71E- 


-02 


1.86987 


Role of JAK family kinases in 


2.52E - 


F 00 


7.14E- 


-02 


1 .79928 


IL-6-type Cytokine Signaling 












Amyloid Processing 


3.20E- 


F 00 


4.92E- 


-02 


1.5744 


Melanocyte Development 


3.73E - 


F 00 


4.21E- 


-02 


1 .57033 


and Pigmentation Signaling 












Gai Signaling 


3.16E - 


F 00 


2.96E- 


-02 


0.93536 


Dopamine Receptor Signaling 


2.70E - 


F 00 


3.12E- 


-02 


0.8424 


CDK5 Signaling 


2.52E- 


F 00 


3.09E- 


-02 


0.77868 


a-Adrenergic Signaling 


2.52E - 


F 00 


2.75E- 


-02 


0.693 


IGF-1 Signaling 


2.40E - 


F 00 


2.80E- 


-02 


0.672 


CNTF Signaling 


1.90E + 00 


3.51E- 


-02 


0.6669 


Neuropathic Pain Signaling In 


2.36E - 


F 00 


2.75E- 


-02 


0.649 


Dorsal Horn Neurons 












G Beta Gamma Signaling 


2.50E- 


F 00 


2.48E- 


-02 


0.62 


IL-2 Signaling 


1.88E -I 


- 00 


3.28E- 


-02 


0.61664 


Thrombopoietin Signaling 


1.85E H 


- 00 


3.12E- 


-02 


0.5772 


Renin-Angiotensin Signaling 


2.26E- 


F 00 


2.38E- 


-02 


0.53788 


Gas Signaling 


2.22E - 


F 00 


2.40E- 


-02 


0.5328 


Role of JAK1 and JAK3 in yc 


1.75E -I 


- 00 


2.94E- 
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0.5145 


Cytokine Signaling 












GM-CSF Signaling 


1.74E -I 


- 00 


2.94E- 


-02 


0.51156 


Clathrin-mediated 


2.49E - 


F 00 


2.02E- 


-02 


0.50298 


Endocytosis Signaling 












Regulation of Cellular 


1.82E H 


- 00 


2.74E- 


-02 


0.49868 


Mechanics by Calpain 












Protease 












CREB Signaling in Neurons 


2.58E- 


F 00 


1.93E- 


-02 


0.49794 


Glutamate Receptor Signaling 


1.79E H 


- 00 


2.78E- 


-02 


0.49762 


Synaptic Long Term 


2.14E + 00 


2.31 E- 


-02 


0.49434 


Potentiation 












Antiproliferative Role of 


1.73E -I 


- 00 


2.78E- 


-02 


0.48094 


Somatostatin Receptor 2 












Agrin Interactions at 


1.68E H 


- 00 


2.86E- 


-02 


0.48048 


Neuromuscular Junction 












JAK/Stat Signaling 


1.69EH 


- 00 


2.82E- 


-02 


0.47658 


Corticotropin Releasing 


2.21 E - 


F 00 


2.07E- 


-02 


0.45747 


Hormone Signaling 












Neurotrophin/TRK Signaling 


1.68E + 00 


2.63E- 


-02 


0.44184 


P2Y Purigenic Receptor 


2.12E - 


F 00 


2.08E- 


-02 


0.44096 


Signaling Pathway 












Renal Cell Carcinoma 


1.65E H 


- 00 


2.53E- 


-02 


0.41745 


Signaling 












Insulin Receptor Signaling 


2.03E + 00 


2.01 E- 


-02 


0.40803 


Melatonin Signaling 


1.65EH 


- 00 


2.47E- 


-02 


0.40755 


FLT3 Signaling in 


1.60E H 


- 00 


2.53E- 


-02 


0.4048 


Hematopoietic Progenitor 












Cells 












GNRH Signaling 


2.05E- 


F 00 


1.96E- 


-02 


0.4018 


cAMP-mediated signaling 


2.23E - 


F 00 


1.77E- 


-02 


0.39471 


Ephrin B Signaling 


1.61 EH 


- 00 


2.44E- 


-02 


0.39284 


Prolactin Signaling 


1.62E H 


- 00 


2.38E- 


-02 


0.38556 


BMP signaling pathway 


1.59EH 


- 00 


2.33E- 


-02 


0.37047 
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Table 6 | Ingenuity Canonical Signaling Pathway analysis of proteins 
differentially regulated in the hippocampus of aged BTBR compared 
to WT controls. 



Canonical pathways -log(p-value) Ratio 10 x Hybrid 



Oxidative phosphorylation 


1.27E - 


-01 


1.25E- 


-01 


15.875 


Mitochondrial dysfunction 


1.33E - 


^01 


8.37E- 


-02 


11.1321 


GABA receptor signaling 


5.19E + 00 


1.07E- 


-01 


5.5533 


Glutamate degradation III (via 


3.02E - 


- 00 


1.67E- 


-01 


5.0434 


4-aminobutyrate) 












GM-CSF signaling 


4.44E - 


- 00 


8.82E- 


-02 


3.91608 


Huntington's disease signaling 


6.46E - 


- 00 


5.16E- 


-02 


3.33336 


G protein signaling mediated 


3.47E - 


- 00 


9.09E- 


-02 


3.15423 


by Tubby 












Melatonin signaling 


4.18E - 


- 00 


7.41 E- 


-02 


3.09738 


RhoA signaling 


4.52E - 


- 00 


6.25E- 


-02 


2.825 


I inirl anfinpn nrpspntatinn hv 


2.68E - 


- 00 


1.00E- 


-01 


2.68 


CD1 












Regulation of actin-based 


3.77E - 


- 00 


6.52E- 


-02 


2.45804 


motility by rho 












Role of NFAT in cardiac 


4.87E - 


- 00 


4.78E- 


-02 


2.32786 


hypertrophy 












CTLA4 Signaling in cytotoxic T 


3.63E - 


h 00 


6.25E- 


-02 


2.26875 


lymphocytes 












TCA Cycle II (Eukaryotic) 


2.90E - 


r 00 


7.50E- 


-02 


2.175 


Breast cancer regulation by 


4.65E - 


h 00 


4.67E- 


-02 


2.17155 


stathminl 












Chemokine signaling 


3.22E - 


h 00 


6.67E- 


-02 


2.14774 


Synaptic long term 


3.72E - 


h 00 


5.38E- 


-02 


2.00136 


potentiation 












Guanine and guanosine 


1.71 E H 


- 00 


1.11 E- 


-01 


1.8981 


salvage 1 












CREB signaling in neurons 


4.21 E - 


h 00 


4.35E- 


-02 


1.83135 


p-alanine degradation 1 


1.71 E H 


- 00 


1.00E- 


-01 


1.71 


Axonal guidance signaling 


5.19E - 


- 00 


3.29E- 


-02 


1.70751 


GNRH signaling 


3.51E - 


- 00 


4.58E- 


-02 


1 .60758 
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2.65E - 


- 00 


5.97E- 


-02 


1 58205 


ARP-WASP complex 












4-Aminobutyrate degradation 1 


1.53E H 


- 00 


1.00E- 


-01 


1.53 


Superpathway of inositol 


3.94E - 


- 00 


3.85E- 


-02 


1.5169 


phosphate compounds 












Creatine-phosphate 


1.32E H 


- 00 


1.11 E- 


-01 


1 .4652 


biosynthesis 












Pyruvate fermentation to 


1.32E H 


- 00 


1.11 E- 


-01 


1 .4652 


lactate 












D-myo-inositol-5-phosphate 


3.32E - 


- 00 


4.32E- 


-02 


1 .43424 


metabolism 












fMLP signaling in neutrophils 


3.13E - 


- 00 


4.55E- 


-02 


1.42415 


Methylglyoxal degradation 1 


1.53E H 


- 00 


9.09E- 


-02 


1.39077 


Clathrin-mediated endocytosis 


3.29E - 


- 00 


4.04E- 


-02 


1.32916 


signaling 












Ephrin receptor signaling 


3.45E - 


- 00 


3.81E- 


-02 


1.31445 


Protein kinase A signaling 


3.95E - 


- 00 


3.18E- 


-02 


1.2561 


Thrombin signaling 


3.20E - 


- 00 


3.79E- 


-02 


1.2128 


3-phosphoinositide 


3.10E H 


- 00 


3.87E- 


-02 


1.1997 


biosynthesis 
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-log(p-value) 


Ratio 


10 x Hybrid 


D-myo-inositol 
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1.19262 
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biosynthesis 
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8.33 E-02 


1.17453 


Neuropathic pain signaling in 


2 50E - 


l oo 


4.59E-02 
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dorsal horn neurons 










G beta gamma signaling 


2.73E - 
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4.13E-02 
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nNOS signaling in neurons 


1 95E -f 


- 00 


5.77 E-02 


1.12515 


Signaling by rho family 


3 24E - 


l oo 


3.36E-02 


1 .08864 


GTPases 










Arginine degradation I 


1 / IC j 

1.4 i c i 


- uu 


7.69E-02 


1 .08429 


(arginase pathway) 










Ephrin B signaling 


2 21 E - 


i_ oo 


4.88E-02 


1 .07848 


Purine nucleotides de novo 


9 OOP 

/.z/t - 


i_ nn 
h UU 


4.76E-02 


1 .05672 


biosynthesis II 










Leptin signaling in obesity 


2 19E - 


- 00 


4.71 E-02 


1.03149 


Semaphorin signaling in 


1 J 
I.OOL -f 


nn 

- UU 


5.56E-02 


1.01748 


neurons 










lnosine-5'-phosphate 


1 53E ^ 


- 00 


6.25E-02 


0.95625 


biosynthesis II 










Actin cytoskeleton signaling 


2.88E - 


i_ oo 


3.31 E-02 


0.95328 


3-phosphoinositide 


z.o / rz - 


i_ nn 
r uu 


3.70E-02 


0.9509 


degradation 










Rac signaling 


2 43 E - 


[_ oo 


3.91 E-02 


0.95013 


iCOS-iCOSL signaling in T 


9 QQC 

z.oot - 


i nn 
h UU 


3.97E-02 


0.94486 


helper cells 










Dopamine receptor signaling 


2 17E - 


- 00 


4.17E-02 


0.90489 


Molybdenum cofactor 


1 /|1C j 

1.4- 1 t - 
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- uu 


6.25E-02 


0.88125 


biosynthesis 










Aldosterone signaling in 


Z.oolz - 
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3. 55 E-02 


0.8449 


epithelial cells 










CD28 signaling in T helper 


2 23E - 


l oo 


3.68E-02 


0.82064 


cells 










B cell receptor signaling 


9 
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i_ nn 
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3.43 E-02 


0.7889 


2-oxobutanoate degradation I 


1.32E -j 


- 00 
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PKC6 signaling in T 


2.23E - 


h 00 


3.47 E-02 


0.77381 


lymphocytes 










Virus entry via endocytic 


1.93E -i 


- 00 


3.96E-02 


0.76428 


pathways 










P2Y purigenic receptor 


2.14E - 


- 00 


3.47 E-02 


0.74258 


signaling pathway 
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0.74 
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biosynthesis 
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0.72225 


in cAMP signaling 










PI3K signaling in B 


2.05E - 


h 00 


3. 50 E-02 


0.7175 


lymphocytes 










a-Adrenergic signaling 


1.95E -i 


- 00 


3.67E-02 


0.71565 


Protein ubiquitination pathway 


2.41 E - 


h 00 


2.96E-02 


0.71336 


Gluconeogenesis 1 


1.67E -f 


- 00 


4.26E-02 


0.71142 
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Canonical pathways -log(p-value) Ratio 10 x Hybrid 



1 1 it ^ m a to rorontnr cmna inn 
UlULdllldLc [fcrL/fcjpLUI blylldllMLj 


1 69E H 


- 00 


4 17E— 02 


0 70473 


U rlULfcJIM UUUpifcJU [fcfUcpLUI 


2 38E - 


f- 00 


2 90E— 02 


0 6902 


signaling 










ILr1 signaling 


1 R7F J 
l.o/ c i 


nn 
- uu 


O.D / f— uz 


n CQR9Q 


Antipro iferative rols of 


I.OUF i 


nn 
- uu 


A 17F n? 

^f. I / F— UZ 


U.DD / Z 


somatostatin rscsptor 2 










Rsmodsling of epithelial 


1.53E H 


- 00 


4.29E-02 


0.65637 


adherens junctions 










JAK/stat signaling 


1.54E H 


- 00 


4.23E-02 


0.65142 


Glioma signaling 


1.84E H 


- 00 


3.54E-02 


0.65136 


Role of NFAT in regulation of 


2.15E - 


f- 00 


3.00E-02 


0.645 


the immune response 










Fatty acid p-oxidation I 


1.45E H 


- 00 


4.44E-02 


0.6438 


RhoGDI signaling 


2.11 E - 


- 00 


2.97E-02 


0.62667 


Cardiac p-adrenergic signaling 


1.96E H 


- 00 


3.16E-02 


0.61936 


Cardiac hypertrophy signaling 


2.15E - 


f- 00 


2.80E-02 


0.602 


Re ax n signs ng 


1.93E H 


- 00 


3.05E-02 


0.58865 


Renal cell carcinoma signaling 


1.50E H 


- 00 


3.80E-02 


0.57 


Calcium signaling 


2.05E - 


h 00 


2.76E-02 


0.5658 


Phospholipase C signaling 


2.06E- 


f- 00 


2.64E-02 


0.54384 


Gotq signaling 


1.80E + 00 


2.92E-02 


0.5256 


Gas signaling 


1.60E + 00 


3.20E-02 


0.512 


Molecular mechanisms of 


2.15E - 


h 00 


2.32E-02 


0.4988 


cancer 










Role of tissue factor in cancer 


1.61 EH 


- 00 


3.08E-02 


0.49588 


Gap junction signaling 


1.71 E H 


- 00 


2.76E-02 


0.47196 


p70S6K signaling 


1.51 eh 


- 00 


3.03E-02 


0.45753 


Androgen signaling 


1.60E H 


- 00 


2.76E-02 


0.4416 


Gai signaling 


1.48E H 


- 00 


2.96E-02 


0.43808 


cAMP-mediated signaling 


1.65E H 


- 00 


2.65E-02 


0.43725 


PPARct/RXRa activation 


1.51 eh 


- 00 


2.50E-02 


0.3775 


Insulin receptor signaling 


1.38E H 


- 00 


2.68E-02 


0.36984 


eNOS signaling 


1.39E H 


- 00 


2.58E-02 


0.35862 


Glucocorticoid receptor 


1.35E H 


- 00 


2.01 E-02 


0.27135 



signaling 

Canonical signaling pathway annotation of the hippocampal proteins possessing 
BTBR:WT iTRAQ expression ratios < 12 or < 0.8 was created using Ingenuity 
Pathway Analysis. Each signaling pathway significantly populated (p < 0.05) by 
at least two independent proteins are depicted. For each populated signaling 
pathway the—log^ of the enrichment probability (p-value), the enrichment ratio 
and the calculated hybrid score [10 * (enrichment ratio * —logio(p-value!)I are 
depicted. 

transmembrane protein 2) and Cep290 (centrosomal protein 
290)] (Figure 6D). 

DISCUSSION 

In the present study we evaluated the behavior, synaptic protein 
alterations and targeted proteomic signatures in the hippocam- 
pus and cortex of 15 month old BTBR and WT mice to determine 
whether the ASD-like phenotype of BTBR mice was resistant to 
neuronal aging. Thus, we assessed both the behavioral phenotype 
and proteomic profiles between aged ASD-like mice (abnormal 
CNS) and aged wild type mice (with a relatively normal CNS). 



We found that several neurosynaptic proteins were significantly 
increased in both the hippocampus and cortex of the aged ASD- 
like mice in comparison to the aged wild type. In particular, 
the phosphorylated form of the pre-synaptic marker synapsin 
1, and the post-synaptic marker spinophilin were increased in 
both the aged BTBR hippocampus and cortex; PSD95 was sig- 
nificantly elevated in aged BTBR hippocampus. Aged BTBR mice 
still presented a behavioral phenotype typically observed in young 
autistic mice, suggesting that the physiological actions underlying 
the ASD-like condition of the BTBR mice is relatively resistant to 
aging. 

The maintenance of the ASD-phenotype into advanced age in 
BTBR mice is intriguing because in the normal developing CNS 
synapses are overproduced but then become pruned in response 
to neuronal activity during post-natal development (Petit, 1988; 
Cohen and Greenberg, 2008; Ebert and Greenberg, 2013). It is 
proposed that autism is caused by dysregulation of the machin- 
ery controlling synaptic messenger RNA translation leading to 
increased synthesis of synaptic proteins and lack of synaptic prun- 
ing, in turn inducing a hyperconnectivity and over-reinforcement 
of neuronal circuits. In fact it is one of the hallmark character- 
istics of ASD observed in the brains of both mouse models and 
humans with ASD (Kelleher and Bear, 2008; Gkogkas et al, 2013). 
Fragile X syndrome, for example which is the leading cause of 
inherited intellectual disability and ASD (Bhakar et al., 2012), is 
caused by mutations in FMR1 gene resulting in decreased expres- 
sion of fragile X mental retardation protein (FMRP). This protein 
modulates the rate of translational elongation of mRNAs that 
encode synaptic proteins (Darnell et al, 2011). A decrease in 
FMRP impairs regulation of synaptic protein expression (Bhakar 
et al, 2012) and post-natal synapse elimination (Harlow et al., 
2010). 

As increased synaptic marker expression is associated with 
the autistic-like mouse phenotype it is not surprising that we 
found that the BTBR T + tf/j mice maintained their reduced 
social behavior into old age (Figure 1). This is consistent with 
the limited human literature that found that social impairments 
in individuals with ASD persist into adulthood (Ballaban-Gil 
et al, 1996; Szatmari, 1998; Howlin and Moss, 2012). Although 
increased synaptic markers have been found to be deleterious 
for the developing CNS, it is possible that the increased synapses 
and hyperconnected neuronal circuits be an effective advantage in 
the aging process where typically nervous system connectivity is 
degraded through accumulated damage. BDNF and TrkB mRNAs 
are localized in dendritic spines, thus age-related spine den- 
sity decreases in older individuals suggest a link between BDNF 
and TrKB mRNA expression and impaired synaptic connections 
with aging (Tapia-Arancibia et al., 2008). Indeed decreases in 
both BDNF and TrkB expression have been documented in the 
hippocampus and in the neocortex during physiological aging 
(Phillips et al, 1991; Murray et al., 1994; Connor et al., 1997; 
Ferrer et al., 1999; Quartu et al., 1999; Hock et al, 2000; Holsinger 
et al, 2000; Fahnestock et al, 2002; Romanczyk et al, 2002; 
Tapia-Arancibia et al, 2008; Zhao et al, 2010). 

We found that levels of mature BDNF were significantly 
decreased in the aged BTBR hippocampus and cortex compared 
with the aged WT (Figures 2 A, 3A), suggesting an accelerated 
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FIGURE 6 | Identification of coherent regulatory protein signatures in 
BTBR CNS tissues. Vennplex diagram analysis of the upregulated (BTBR:WT 
iTRAQ ratio > 1.2) or downregulated (iTRAQ ratio < 0.8) proteins in both the 
hippocampus and cortex are depicted (A). A small coherently-regulated 
protein subset across both tissues was evident and comprised a focused 
protein set strongly linked to human ASD-like conditions (B). 7exrraus/-based 



collective processing interrogation of the coherently-regulated BTBR protein 
set revealed, in the resultant hierarchical word-cloud, a strong potential role 
of skeletal modeling activity in this protein set (C). This skeletal signaling 
dependence of ASD was confirmed using the individual processing module 
of Textrous! (D). In this heatmap output correlation strength between protein 
and functional term is indicated by the intensity of the teal-colored blocks. 



aging phenomena, but did not find a significant change in the 
overall level of active TrkB receptor (Figures 2D, 3D), in the 
BTBR mice. In addition, we also found that there was a non- 
significant trend for the increase in the expression in both cortical 
and hippocampal levels of Aktl, one of the most important 
downstream targets of TrkB signaling. Thus, via a potentially 
complex compensatory mechanism, the BTBR mice may have 
increased their BDNF signaling efficiency and downstream signal 
transduction capacity compared to WT mice. 

Our findings of reduced BDNF in both the hippocampus and 
cortex in old BTBR mice compared to old WT mice despite seem- 
ingly intact memory functioning may be related to this efficiency 
potentiation of the BDNF-TrkB system in the BTBR mice. These 
effects may be related to changes in endocytic proteins such as 
Itsnl, Golga2 and Ap3b2 (Table 7) that may affect the signaling 
dynamics of the TrkB receptor (Mattson et al., 2004). In addition 



it is possible that the elevated levels of various synaptic mark- 
ers served to protect the aging BTBR brain against age-related 
declines in cognition often associated with decreased levels of 
BDNF. 

Given the complexity of the ASD, it is unlikely that BDNF 
signaling and various neuronal cell markers are the only pro- 
tein products altered in the presence of ASD; thus we performed 
global quantitative iTRAQ proteomics, coupled to advanced 
bioinformatics in order to gain a comprehensive, systems-level 
understanding of potential changes in aged BTBR mice compared 
to WT counterparts. With respect to current relevant literature we 
discovered that many of these proteins differentially regulated in 
the brain of the old BTBR mice demonstrated strong correlations 
with those altered in humans with neurodevelopmental disabili- 
ties (Table S7). Nearly twenty percent of the coherently-regulated 
proteins across the hippocampus and cortex of aged BTBR mice 
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also are strongly linked to human ASD phenotypes, suggest- 
ing that future investigation of the remaining factors may prove 
extremely fruitful for autism-related research (Figure 6B). With 
respect to the highly-focused, coherently-regulated protein set 
(Figure 6B: Table 7) Prrt2 was coherently altered in both the hip- 
pocampus and cortex of aged BTBR mice compared to aged WT 
(Table 7). Alterations in activity of Prrt2 have been found in the 
frontal cortex of individuals with autism (Ji et al., 2012) and have 
also been genetically-associated with ASD (Weber et al., 2013). 
Vatll was up-regulated in both the hippocampus and cortex of 
the old BTBR mice. 

We performed functional informatics data clustering to gain 
a wider appreciation of the signaling phenotype present in aged 
BTBR mice. Using GO term annotation it was evident in the 
cortex that considerable alterations in ER (endoplasmic reticu- 
lum), vesicular transport, cytoskeletal remodeling, and energy 
generation were apparent (Figures 4E,F). The functional GO 
term clustering of differentially expressed proteins in the BTBR 
hippocampus also indicated a role of energy generation, neu- 
rotransmitter release and synaptic activity but did not display 
a robust skeletal remodeling component (Table 2). We also 
employed KEGG and Ingenuity Canonical Signaling pathway 
analysis for the differentially-regulated cortex and hippocam- 
pal protein sets. At this higher level of functional interrogation 
(compared to GO term clustering) we also observed a subtle 
divergence in tissue behavior, i.e., cortical KEGG pathway anal- 
ysis demonstrated a strong bias toward energy metabolic and 
neurotrophic/neuroprotective mechanisms (Figure 5A) while the 
hippocampus demonstrated a more degenerative phenotype 
associated with metabolic instabilities and calcium management 
factors (Figure 5B). Interestingly and to highlight the need for 
multidimensional informatics analyses for physiological interpre- 
tation, the canonical signaling pathway analysis of both tissue 
datasets revealed additional processes that could affect the ASD 
phenotype, i.e., the prominence of cytokine signaling activity in 
the cortex and the identification of the potential role of melatonin 
function in both CNS tissues (Figures 5C,D). 

To further reinforce the importance of applying unbiased 
informatics for the analysis of complex data associated with a 
highly nuanced phenotype such as ASD we found that analysis 
of different subgroups of the data yielded interesting, tissue- 
specific insights into the effects of aging on ASD. Using our 
novel Textrous! platform we were able to identify differential 
tissue activity, i.e., enhanced glutamatergic cortical activity coin- 
cident with attenuated kinase signaling (Figure SI). These cor- 
tical changes occurred contemporaneously with an elevation of 
hippocampal protein misfolding/aggregation coincident with an 
attenuated CNS myelination and oligodendrocyte functional- 
ity (Figure S2). Both of these hippocampal processes have been 
recently associated with dysfunctional brain development (Mitew 
etal., 2013). 

When approaching complex issues such as ASD with a 
systems-biology approach it is important to maintain this concept 
of multiple levels of protein interaction, function and hierarchi- 
cal architecture, i.e., functions can be tissue specific, expression 
polarity specific while also existing contemporaneously at a level 
above tissue-or polarity-dependence. For example we found that 



Table 7 | VennPlex Venn diagram analysis of up and downregulated 
proteins across the cortex and hippocampus of BTBR mice compared 
to WT controls. 



Protein 


Cortex up 


Cortex down 


Hippo up 


Hippo down 


PRRT2 


1.951168338 




2.381953722 




DYNLL1 


1.672920918 




1.737545789 




DBNL 


1.626693344 




1.24548243 




AP3B2 


1.581785753 




1.789312776 




VAT1L 


1.391942674 




1.401775589 




WDR33 


1.379833997 




2.510567518 




GOLGA2 


1.292500223 




1.60040989 




LANCL1 


1.279121599 




1.320074872 




ME3 


1.270939623 




1.62730561 




ITSN1 


1.269638936 




1.774926636 




CEP290 


1.238621608 




1.454207016 




PCYOX1 L 


1.209317219 




1.860882083 




ERMN 




-1.243267187 




-1.500264003 


CSRP1 




-1.365689412 




-1.509867503 


DB 




-1.496288514 




-2.056916039 


NUCB1 




-1.631205729 




-3.561025937 


EIF5A 


1.517756095 






-1.65673365 


MMAA 


1.391228352 






-1.887879989 


PTPN11 


1.247745108 






-3.093767449 


ZXDB 


1.214676409 






-1.393599577 


SLC17A7 


1.204082729 






-1.315805498 


DPYSL4 


1.195193712 






-1.314442062 


OAT 




-1.307905257 


1.362629072 




Ptcd 


1.828557315 








PCBP2 


1.792502849 








QDPR 


1.7625157 








SRSF4 


1.533013648 








ADCY2 


1.511246891 








NDUFS4 


1.47998035 








SLC9A3R1 


1.473396618 








LETM1 


1.455678341 








PRDX3 


1.44243043 








IDH2 


1.441010435 








EXOC5 


1.41342868 








HNRNPUL2 


1.405691244 








PTCD1 


1.371753507 








DDX5 


1.365841661 








PPP1R9B 


1.346708939 








DPP6 


1.340020376 








LAP3 


1.303333413 








KCNAB2 


1.296085587 








SLC12A5 


1.252190507 








HOOK3 


1.250023514 








COX5A 


1.247002966 








WDR1 


1.244912709 








GRM2 


1.243617897 








CAB39 


1.2340405 








RTN3 


1.233971862 








TF 


1.233243832 








FBX041 


1.231290747 








RAB7A 


1.226744118 








RTN4 


1.216638369 








TCP1 


1.215498024 
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Protein Cortex up Cortex down Hippo up Hippo down 



PRKAR2B 


1.215235782 






CCT3 


1.20748828 






PRDM5 


1.198502382 






COQ9 




-1.24961385 




CLDN11 




-1.260614004 




NCALD 




-1.266374541 




F5 




-1.269955566 




EXOC6B 




-1.275380087 




CTTN 




-1.287911813 




VBP1 




-1.361766656 




PICALM 




-1.466374966 




CAPNS1 




-1.756445493 




TOMM20 




-1.819905488 




MAPK3 




-1.918813356 




NUDC 






1.732009523 


RAP2B 






1.726355193 


KIF2A 






1.663032474 


SETDB1 






1.648556725 


ATP6V1 D 






1.640964877 


SYNGR1 






1.58532711 


DPYSL3 






1.582802831 


NCAN 






1.571585688 


CTNND2 






1.562506309 


CNRIP1 






1.554254382 


UTP20 






1.544855143 


ME1 






1.478800627 


PFN1 






1.477145399 


PPME1 






1.455765369 


APP 






1.455465641 


LANCL2 






1.455132031 


GNAQ 






1.453050806 


ARPC5L 






1.451224208 


PRDX2 






1.448929334 


CYFIP1 






1.432347642 


HSPA2 






1.431503207 


SFPQ 






1.427372534 


TXLNG 






1.42119969 


GRB2 






1.4157442 


CPNE6 






1.414165286 


PDXP 






1.410098436 


GNB2 






1.408516442 


GLOD4 






1.391614359 


OXR1 






1.380069445 


VAMP2 






1.373709575 


HSPA1A 






1.373028035 


HSPA1B 






1.373028035 


PSD3 






1.371411227 


BSG 






1.369172074 


PFN2 






1.368469072 


ATP5J2 






1.368321641 


GNB1 






1.366417736 


DLG4 






1.363981591 


SH3GLB2 






1.36363306 


KTN1 






1.362629072 


RAB3A 






1.362173026 



(Continued) 
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Protein Cortex up Cortex down 


Hippo up Hippo down 


CCT7 


1.3567651 


SDHC 


1.347654394 


BCAN 


1.342898167 


HSPH1 


1.337373319 


ADD2 


1.336930818 


AP2A2 


1.335723865 


DYNLL2 


1.322647505 


GDA 


1.312825433 


TKT 


1.3105561 


COX7A2 


1.308575621 


SNCA 


1.305570043 


BAIAP2 


1.303596608 


SH3BGRL3 


1.302583328 


NPM1 


1.301978732 


NIPSNAP1 


1.299090793 


ARPC4 


1.296939524 


NAPA 


1.292138971 


CALB1 


1.290893833 


CAMKV 


1.286413462 


COX6C 


1.286388519 


COX5B 


1.285714649 


NDUFA9 


1.285358088 


ATP6V0D1 


1.275603048 


NSFL1C 


1.274526051 


LDHA 


1.269135901 


ATP6V1 F 


1.266693513 


CLTB 


1.265664013 


SRSF4 


1.261593119 


LRGUK 


1.260891479 


NDUFS8 


1.260447912 


PPP2CB 


1.259397109 


COX4I1 


1.256849924 


CAMK2D 


1.256848229 


DNAJC6 


1.255985647 


COL6A2 


1.253698992 


CRYM 


1.251866833 


PPP3CA 


1.249125255 


SLC2A1 


1.24839376 


PPFIA3 


1.244596343 


CAMK2A 


1.241519397 


SYP 


1.237859198 


AP2B1 


1.237413188 


USP24 


1.233820239 


FABP5 


1.233314773 


SEPT9 


1.226916013 


EHD1 


1.226344591 


CAP2 


1.223807885 


SLC25A3 


1.221941931 


L.MIVI IxZD 


i .zzuoo^+yy 


PTGES3 


1.216665914 


PRDX1 


1.215985447 


ATP5I 


1.215521395 


ATIC 


1.213849668 


GSTM5 


1.213704207 


NDUFS7 


1.211966939 



(Continued) 
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Protein Cortex up 


Cortex down Hippo up 


Hippo down 


DCLK1 


1.207992461 




GPD1L 


1.207286149 




HNRNPL 


1.205847041 




NAPB 


1.202833155 




RTN4 


1.201798491 




PCSK1N 




-1.245474191 


WDR7 




-1.245967263 


MAG 




-1.248958157 


DDB1 




-1.255013637 


MOBP 




-1.256321905 


GAD2 




-1.26014124 


NA 




-1.265874043 


TPT1 




-1.268669401 


IDH3G 




-1.270175029 


PCCB 




-1.272883353 


MAPRE2 




-1.274095945 


MRPS36 




-1.275342647 


MARCKS 




-1.280485807 


PSMC6 




-1.281124693 


ABAT 




-1.283265073 


NTM 




-1.284285221 


CA2 




-1.284359839 


HPCAL1 




-1.284463003 


MOG 




-1.295178624 


PSMA6 




-1.297701354 


SDHA 




-1.300775393 


COPB1 




-1.300938943 


NONO 




-1.301095086 


TUFM 




-1.311237761 


ADSS 




-1.311542644 


IGSF8 




-1.312784192 


COX5A 




-1.313555144 


UQCRFS1 




-1.316410553 


CKB 




-1.318123142 


PIP4K2B 




-1.324477494 


NEFM 




-1.330302859 


PSMD6 




-1.330311752 


NEFL 




-1.344530231 


NDUFA8 




-1.348309239 


BIN1 




-1.354260702 


PTMS 




-1.379162453 


PHB2 




-1.381012946 


NDUFS2 




-1.382654847 


DSTN 




-1.396469687 


CCT5 




-1.396711507 


S100B 




-1.403153034 


CLASP2 




-1.403192685 


EPB41L1 




-1.419486383 


rrl\L 




1 A OA ORQOQ/l 


PSAP 




-1.429458675 


CPLX1 




-1.437904339 


SERPINA3K 




-1.462482626 


AHSA1 




-1.482509362 


SET 




-1.494218612 


ATL1 




-1.49652324 
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Protein Cortex up Cortex down 


Hippo up Hippo down 


PRKACB 


-1.501341659 


NME2 


-1.502355851 


ITPKA 


-1.515605888 


NCS1 


-1.519434163 


EEF1A1 


-1.528884517 


PLCB1 


-1.530645112 


ENPP6 


-1.542384655 


SYNE1 


-1.544411246 


PCP4 


-1.554348855 


MBP 


-1.583600552 


NCEH1 


-1.612308156 


ACSL6 


-1.612342183 


GL01 


-1.659621805 


PPP1R1B 


-1.668571078 


HSD17B10 


-1.719780931 


SIRT2 


-1.776727737 


br nlM 


— I . o I UZ4oyU I 


MORC1 


-1.853643692 


SLC6A11 


-1.877132951 


PDHX 


-1.878732243 


PTK6 


-1.89800334 


DPYSL5 


-1.903465278 


HPRT1 


-2.333095423 



For each specific protein studied fold change values, indicated in the table 
below, were generated from the iTRAQ expression ratios of (> 12 or <0.8) pro- 
teins differentially regulated in the BTBR tissues [Cortex, Hippocampus (Hippo)l 
compared to WT. 



inspection of the coherently-regulated multi-tissue protein set 
revealed an exciting and potentially novel functional association 
of neuronal Wasp (Wiskott-Aldrich syndrome protein) activity 
and ASD (Figures 6C,D). To our current knowledge this is the 
first example of such a functional correlation, however it is rela- 
tively simple to predict that the molecular activities of this skeletal 
remodeling factor could impact a disorder such as ASD in which 
a less synaptically plastic neurological system is present. 

In conclusion, we found that both the social behavioral phe- 
notype and the CNS proteomic profile of this ASD model are 
relatively age-resistant. Aged BTBR mice possessed elevations in 
hippocampal and cortical synaptic markers that likely play a 
role in the etiology of the autistic-like phenotype. Low levels 
of BDNF have been found inconsistently in the ASD literature 
but more consistently in the aging and Alzheimer Disease lit- 
erature (Chadwick et al., 2010b). Diminished mature BDNF in 
the aged BTBR hippocampus/cortex were seen along with com- 
plementary alterations in the downstream signaling pathways 
compared to WT aged mice. These findings in-part reinforce the 
potential long-term neuroprotective capacity of the ASD pheno- 
type with respect to cognitive aging. It is therefore interesting to 
posit that despite lower levels of neurosynaptic plasticity that the 
more robust and reinforced neuronal network in ASD pheno- 
types may indeed act prophylactically over long periods of time 
compared to normal CNS networks that degrade due to unre- 
paired damage and synaptic loss. The complex proteotype of the 
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Table 8 | Textrousl-based collective analysis of coherently-regulated 
BTBR-specific proteins. 



Word 


Cosine similarity 


Z-score 


p-Value 


Motor 


0.524476314 


1.869598279 


0 


.030741909 


Nonsense 


0.478414897 


1.693752796 


0 


.045132642 


Atrophies 


0.4700895 


1.661969501 


0 


.048256387 


Gtpases 


0.468206051 


1.65477919 


0 


048962249 


Exons 


0.466743056 


1.649194013 


0 


.049573817 


Growth-associated 


0.455961506 


1.608034038 


0 


.053917589 


Inherited 


0.455845768 


1.607592194 


0 


.053917589 


Gtpase 


0.455708175 


1.607066915 


0 


054027184 


Autosomal 


0.453900108 


1.600164382 


0 


054799292 


Polymerization 


0.453672232 


1.599294435 


0 


054910301 


Optic 


0.444664336 


1.564905613 


0 


.058791455 


Exon 


0.437675556 


1.538225034 


0 


.062024307 


Gtp-binding 


0.435679912 


1.5306064 


0 


062884696 


Maps 


0.427384101 


1.498936057 


0 


.066936816 


Bundles 


0.426675363 


1.496230358 


0 


.067326828 


Disabilities 


0.423470554 


1.483995579 


0 


06890446 


Gdi 


0.419726975 


1.469703976 


0 


.070780877 


Brothers 


0.419534775 


1.46897023 


0 


.070916395 


Movement 


0.417518799 


1.461273977 


0 


.072007721 


Acuity 


0.417184158 


1.45999644 


0 


.072145037 


Cytoskeleton 


0.414147875 


1.448405032 


0 


.073808525 


n-wasp 


0.412898356 


1.443634832 


0 


074369488 


Frameshift 


0.411123648 


1.436859653 


0 


075358997 


Recessive 


0.409375594 


1.430186227 


0 


07635851 


Wiskott-aldrich 


0.407779814 


1.424094126 


0 


.077223236 


Actin 


0.406644694 


1.419760659 


0 


.077803841 


Wasp 


0.405802073 


1.416543842 


0 


.078241464 


Gtpase-activating 


0.402282931 


1.403109057 


0 


.080308419 


Recycling 


0.395770067 


1.378245348 


0 


084101 644 


Reorganization 


0.391991869 


1.363821581 


0 


.086283783 


Prenylation 


0.389968302 


1.356096352 


0 


.087549584 


Splice 


0.388453273 


1.350312529 


0 


.088507991 


Missense 


0.388195935 


1.349330109 


0 


.088668483 


Streaming 


0.387806369 


1.347842888 


0 


.088829191 


Apparatus 


0.3874916 


1.346641219 


0 


.088990116 


Pericentriolar 


0.387010855 


1.344805911 


0 


.089312617 


Families 


0.386042261 


1.341108175 


0 


.089960226 


Plexiform 


0.385574635 


1.339322953 


0 


.090285336 


Polysomes 


0.383798601 


1.332542712 


0 


091265902 


Guanosine 


0.383107191 


1.329903161 


0 


091759136 


Splicing 


0.382666221 


1.328219703 


0 


.092089053 


Cytoskeletal 


0.381993528 


1.325651608 


0 


.092419849 


Mutational 


0.381893517 


1.325269803 


0 


092585576 


Coiled-coil 


0.381550118 


1.323958833 


0 


.092751522 


gaps 


0.380853905 


1.321300948 


0 


.093250682 


Depolymerizing 


0.380444571 


1.319738261 


0 


.093417509 


1 ntellectual 




I.O I ZZoOZUz. 


0 


094760067 


Mutations 


0.378323187 


1.311639602 


0 


.094760067 


Protrusions 


0.377199886 


1.307351254 


0 


.095606356 


Neurofibromin 


0.376699229 


1.305439931 


0 


.095946424 


Actin-binding 


0.375976678 


1.302681498 


0 


.096287381 


Atrophy 


0.375377983 


1.3003959 


0 


.096800485 


Dynamics 


0.373862886 


1.29461182 


0 


097660115 


Pigmentary 


0.373853033 


1.294574204 


0 


.097660115 



(Continued) 



Table 8 | Continued 



Word 


Cosine similarity 


Z-score 


p-Value 


Neuron 


0.370913941 


1.283353835 


0.099746038 


Multidisciplinary 


0.369651835 


1.27853558 


0.100448528 


Destination 


0.369566192 


1.278208627 


0.100624713 


Cilia 


0.368620901 


1.274599856 


0.101154621 


Gtp-dependent 


0.367048826 


1.268598253 


0.102220532 


Phosphatidic 


0.366734218 


1.267397195 


0.102577645 


Identification 


0.366608729 


1.266918125 


0.102577645 


Disability 


0.365764506 


1.263695196 


0.103115013 


Muscular 


0.365741802 


1.26360852 


0.103115013 


Ciliary 


0.364952202 


1.260594116 


0.103654424 


Meshwork 


0.364710238 


1.259670389 


0.103834681 


Partners 


0.364390606 


1.258450149 


0.104195878 


Inheritance 


0.364260066 


1.257951798 


0.104195878 


Scanning 


0.363734162 


1.255944089 


0.104557985 


f-actin 


0.363367304 


1.254543564 


0.10473938 


Beta-galactoside 


0.36316993 


1.253790062 


0.104921003 


Anonymous 


0.361144496 


1.2460577 


0.106382198 


Protrusion 


0.359889734 


1.241267484 


0.107302874 


Sialyltransferase 


0.359831442 


1.241044947 


0.107302874 


Neurofibromatosis 


0.359791096 


1.24089092 


0.107302874 


Lipofuscin 


0.359746522 


1.240720753 


0.107302874 


Radixin 


0.35878958 


1.237067503 


0.108043541 


Proteins 


0.358631641 


1.236464548 


0.108229282 


Squid 


0.358474059 


1.23586296 


0.108229282 


Screen 


0.356102317 


1.226808525 


0.109911295 


Filopodia 


0.355914644 


1.226092057 


0.110099338 


Extensions 


0.355265349 


1.223613288 


0.110476114 


Lowe 


0.355018906 


1.222672462 


0.110664848 


Pre-mrna 


0.354335273 


1.220062601 


0.111232437 


Leber 


0.354284348 


1.219868188 


0.111232437 


Counseling 


0.354132651 


1.219289065 


0.111422096 


x-l inked 


0.354046801 


1.218961323 


0.111422096 


Cones 


0.353920096 


1.218477608 


0.111611986 


Contig 


0.353576868 


1.217167292 


0.111802108 


Cryomicroscopy 


0.35310173 


1.215353389 


0.112183046 


Amaurosis 


0.353046083 


1.215140949 


0.112183046 


Methanothermobacter 


0.352545212 


1.21322881 


0.11256491 


Neurofibromas 


0.35222224 


1.211995824 


0.11275619 


Endocytic 


0.352220725 


1.211990038 


0.11275619 


Pumping 


0.351073392 


1.207609946 


0.113523631 


Tail 


0.350586481 


1.205751097 


0.113908745 


Dynamin 


0.350367829 


1.204916365 


0.114101651 


Syndrome 


0.349521708 


1.20168619 


0.114681764 


Endosomes 


0.347991562 


1.195844656 


0.11584828 


Blindness 


0.34724489 


1.19299414 


0.116434687 


Visual 


0.347190286 


1.192785679 


0.116434687 



Cosine similarity scores, Z-scores and probability values (p-Value) were calcu- 
lated using collective processing of the 16 coherently-regulated BTBR-specific 
proteins in the cortex and hippocampal datasets. 

BTBR mice revealed molecular and functional signatures highly 
reminiscent of those found in humans with neurodevelopmen- 
tal disabilities including ASD. Many of these identified processes 
play a role in neurodevelopment based on studies in murine mod- 
els, and give more evidence for proteins that are playing a role 
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in neuroprotection allowing for maintenance of cognition in the 
aged BTBR mice. Therefore, it is possible that the elevated levels 
of various synaptic markers and proteins identified proteomi- 
cally such as Picalm, Itsnl, Mobp, Ina, as well as signaling actions 
involving Wasp functionality served to protect the aging BTBR 
brain against age-related declines in cognition often associated 
with decreased BDNE In the current study we observed a robust 
maintenance of cognitive ability despite advanced age that may 
be due to unique neurosynaptic architecture and activity in the 
ASD phenotype. The role of the associated proteins and synaptic 
markers provide a foundation for further investigation into this 
intriguing mechanism. 
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Figure S1 | Collective Textrous! analysis of differentially-regulated cortical 
protein subsets. Functional hierarchical word-clouds were generated from 
the upregulated (BTBR:WT iTRAQ ratio > 1.2) (A) and downregulated 
(BTBR:WT iTRAQ ratio < 0.8) cortical protein data sets. 

Figure S2 | Collective Textrous! analysis of differentially-regulated 
hippocampal protein subsets. Functional hierarchical word-clouds were 
generated from the upregulated (BTBR:WT iTRAQ ratio >1.2) (A) and 
downregulated (BTBR:WT iTRAQ ratio < 0.8) hippocampal protein data 
sets. 

REFERENCES 

American Psychiatric Association. (2012). Diagnostic and Statistical Manual of 
Mental Disorders: DSM-V. Washington, DC: American Psychiatric Association. 

Anitha, A., Nakamura, K., Thanseem, I., Yamada, K., Iwayama, Y., Toyota, 
T., et al. (2012). Brain region- specific altered expression and association of 
mitochondria-related genes in autism. Mol. Autism 3:12. doi: 10.1186/2040- 
2392-3-12 

Ballaban-Gil, K., Rapin, L., Tuchman, T., and Shinnar, S. (1996). Longitudinal 
examination of the behavioral, language, and social changes in a population of 
adolescent and young adults with autistic disorder. Pediatr. Neurol. 15, 217-223. 
doi: 10.1016/S0887-8994(96)00219-6 

Bhakar, A. L., Dolen, G., and Bear, M. F. (2012). The pathophysiology of fragile X 
(and what it teaches us about synapses). Ann. Rev. Neurosci. 35, 417-443. doi: 
10.1146/annurev-neuro-060909-153138 

Blednov, Y A., Walker, D. L., Iyer, S. V., Homanics, G., and Harris, A. R. (2010). 
Mice lacking Gad2 show altered behavioral effects of ethanol, flurazepam and 
gabaxadol. Addict. Biol. 15, 45-61. doi: 10.1111/j.l369-1600.2009.00186.x 

Bohlen, M., Cameron, A., Metten, P., Crabbe, J. C, and Wahlsten, D. (2009). 
Calibration of rotational acceleration for the rotarod test of rodent motor coor- 
dination./. Neurosci. Methods 178, 10-14. doi: 10.1016/j.jneumeth.2008. 11.001 

Bolivar, V. J., Walter, S. R., and Phoenix, J. L. (2007). Assessing autism-like behavior 
in mice: variations in social interactions among inbred strains. Behav. Brain Res. 
176:2126. doi: 10.1016/j.bbr.2006.09.007 

Cai, H., Chen, H., Yi, T, Daimon, C, Boyle, J., Peers, C, et al. (2013). VennPlex- 
a novel Venn diagram program for comparing and visualizing datasets with 



differentially regulated datapoints. PLoS ONE 8:e53388. doi: 10.1371/jour- 
nal.pone.0053388 

Centers for disease control and prevention. (2010).Available online at: http://www. 
cdc.gov/ncbddd/dd/default.html [Accessed, 2013] 

Chadwick, W., Brenneman, R., Martin, B., and Maudsley, S. (2010a). Complex 
and multidimensional lipid raft alterations in a murine model of Alzheimer's 
disease. Int. }. Alzheimer's Disease 2010:604792. doi: 10.4061/2010/604792 

Chadwick, W., Mitchell, N., Caroll, J., Zhou, Y., Park, S. S., Wang, L., et al. (2011). 
Amitriptylinemediated cognitive enhancement in aged 3 x Tg Alzheimer's dis- 
ease mice is associated with neurogenesis and neurotrophic activity. PLoS ONE 
6:e21660. doi: 10.1371/journal.pone.0021660 

Chadwick, W., Zhou, Y, Park, S. S., Wang, L., Mitchell, N., Stone, M., et al. (2010b). 
Minimal peroxide exposure of neuronal cells induces multifaceted adaptive 
responses. PLoS ONE 5:el4352. doi: 10.1371/journal.pone.0014352 

Chao, H. T, Chen, H., Samaco, R. C, Xue, M., Chahrour, M., Yoo, J., et al. 
(2010). Dysfunction in GABA signaling mediates autism-like stereotypies and 
Rett syndrome phenotypes. Nature 468, 263-269. doi: 10.1038/nature09582 

Chen, H., Martin, B., Daimon, C, Siddiqui, S., Luttrell, L., and Maudsley, S. 

(2013) . Textrous!: extracting semantic textual meaning from gene sets. PLoS 
ONE 8:e62665. doi: 10.1371/journal.pone.0062665 

Cohen, S., and Greenberg, M. E. (2008). Communication between the synapse and 
the nucleus in neuronal development, plasticity, and disease. Annu. Rev. Cell 
Dev. Biol. 24, 183-209. doi: 10.1146/annurev.cellbio.24.110707.175235 

Connor, B., Young, D., Yan, Q., Faull, R. L., Synek, B., and Dragunow, M. (1997). 
Brain-derived neurotrophic factor is reduced in Alzheimer's disease. Mol. Brain 
Res. 49:7181. doi: 10.1016/S0169-328X(97)00125-3 

Coppieters, F., Casteels, I., Meire, R, De Jaegere, S., Hooghe, S., van Regemorter, N., 
et al. (2010). Genetic screening of LCA in Belgium: predominance of CEP290 
and identification of potential modifier alleles in AHI1 of CEP290-related 
phenotypes. Hum. Mutat. 31, E1709-E1766. doi: 10.1002/humu.21336 

Crawley, J. N. (2007). Mouse behavioral assays relevant to the symptoms of autism. 
BrainPathol. 17, 448-459. doi: 10.1111/j.l750-3639.2007.00096.x 

Cukier, H. N., Dueker, N., Slifer, S., Lee, J. M., Whitehead, P. L., Lalanne, E., et al. 

(2014) . Exome sequencing of extended families with autism reveals genes shared 
across neurodevelopmental and neuropsychiatric disorders. Mol. Autism 5:1. 
doi: 10.1186/2040-2392-5-1 

Cukier, H. N., Lee, J. M., Ma, D., Young, J. I., Mayo, V., Butler, B. L., et al. (2012). 
The expanding role of MBD genes in autism: identification of a MECP2 dupli- 
cation and novel alterations in MBD5, MBD6, and SETDB1. Autism Res. 5, 
385-397. doi: 10.1002/aur.l251 

Darnell, J. C, Van Driesche, S. J., Zhang, C, Hung, K. Y, Mele, A., Fraser, C. E., 
et al. (2011). FMRP stalls ribosomal translocation on mRNAs linked to synaptic 
function and autism. Cell 146, 247-261. doi: 10.1016/j.cell.2011.06.013 

Defensor, E. B., Pearson, B. L., Pobbe, R. L. H., Bolivar, V. J., Blanchard, D. C., and 
Blanchard, R. J. (2011). A novel social proximity test suggests patterns of social 
avoidance and gaze aversion-like behavior in BTBR T + tf/J mice. Behav. Brain 
Res. 217, 302-308. doi: 10.1016/j.bbr.2010.10.033 

Dere, E., Huston, J. P., and De Souza Silva, M. A. (2007). The pharma- 
cology, neuroanatomy and neurogenetics of one-trial object recognition in 
rodents. Neurosci. Biobehav. Rev. 31, 673-704. doi: 10.1016/j.neubiorev.2007. 
01.005 

Dickson, T. C, Chuckowree, J. A., Chuah, M. I., West, A. K., and Vickers, 
J. C. (2005). Alphainternexin immunoreactivity reflects variable neuronal 
vulnerability in Alzheimer's disease and supports the role of the beta- 
amyloid plaques in inducing neuronal injury. Neurohiol. Dis. 18, 286-295. doi: 
10.1016/j.nbd.2004.10.001 

Ebert, D., and Greenberg, M. (2013). Activity- dependent neuronal signaling and 
autism spectrum disorder. Nature 493, 327-337. doi: 10.1038/naturell860 

Fahnestock, M., Garzon, D., Holsinger, R. M., and Michalski, B. (2002). 
Neurotrophic factors and Alzheimer's disease: are we focusing on the wrong 
molecule? /. Neural. Transm. Suppl. 62, 241-252. doi: 10.1007/978-3-7091-6139- 
5_22 

Ferrer, I., Marin, C, Rey, M. J., Ribalta, T, Goutan, E., Blanco, R., et al. (1999). 
BDNF and full-length and truncated TrkB expression in Alzheimer disease: 
implications in therapeutic strategies. /. Neuropathol. Exp. Neurol. 58, 729-739. 
doi: 10.1097/00005072-199907000-00007 

Feyder, M., Karlsson, R., Mathur, P., Lyman, M., Bock, R., Momenan, R., et al. 
(2010). Association of mouse Dlg4 (PSD-95) gene deletion and human DLG4 
gene variation with phenotypes relevant to autism spectrum disorders and 



Frontiers in Aging Neuroscience 



www.frontiersin.org 



September 2014 | Volume 6 | Article 225 | 23 



Jasien et al. 



Aging and autism spectrum disorder 



Williams' syndrome. Am. J. Psychiatry 167, 1508-1517. doi: 10.1176/appi.ajp. 
2010.10040484 

Fombonne, E., Simmons, H., Ford, T., Meltzer, H., and Goodman, R. (2001). 
Prevalence of pervasive developmental disorders in the British nationwide sur- 
vey of child mental health. /. Am. Acad. Child Adolesc. Psychiatry 40, 820-827. 
doi: 10. 1097/00004583-200 107000-000 1 7 

Fruhmesser, A., Blake, J., Haberlandt, E., Baying, B., Raeder, B., Runz, H., et al. 
(2013). Disruption of EXOC6B in a patient with developmental delay, epilepsy, 
and a de novo balanced t(2;8) translocation. Eur. J. Hum. Genet. 21, 1 177-1 180. 
doi: 10.1038/ejhg.2013.18 

Gillberg, C, and Wing, L. (1999). Autism: not an extremely rare disorder. Acta 
Psychiatr. Scand. 99, 399-406. doi: 10.1111/j.l600-0447.1999.tb00984.x 

Gkogkas, C. G., Khoutorsky, A., and Ran, I. (2013). Autism-related deficits via 
dysregulated eIF4E-dependent translational control. Nature 493, 371-377. doi: 
10.1038/naturell628 

Griswold, A. J., Ma, D., Cukier, H. N., Nations, L. D., Schmidt, M. A., Chung, R. H., 
et al. (2012). Evaluation of copy number variations reveals novel candidate genes 
in autism spectrum disorder-associated pathways. Hum. Mol. Genet. 21:3523. 
doi: 10.1093/hmg/ddsl64 

Harlow, E. G., Till, S. M., Russell, T. A., Wijetunge, L. S., Kind, P., and Contractor, 
A. (2010). Critical period plasticity is disrupted in the barrel cortex of FMR1 
knockout mice. Neuron 65, 385-398. doi: 10.1016/j.neuron.2010.01.024 

Harold, D., Abraham, R., Hollingworth, P., Sims, R., Gerrish, A., Hamshere, M. 
L., et al. (2013). Genome-wide association study identifies variants at CLU and 
PICALM associated with Alzheimer's disease. Nat. Genet. 45, 712. doi: 10.1038/ 
ng0613-712a 

Hedges, D. J., Hamilton-Nelson, K. L., Sacharow, S. J., Nations, L., Beecham, G. 
W, Kozhekbaeva, Z. M., et al. (2012). Evidence of novel fine-scale structural 
variation at autism spectrum disorder candidate loci. Mol. Autism 3:2. doi: 
10.1186/20402392-3-2 

Hock, C, Heese, K., Hulette, C, Rosenberg, C, and Otten, U. (2000). Region- 
specific neurotrophin imbalances in Alzheimer disease: decreased levels of brain 
derived neurotrophic factor and increased levels of nerve growth factor in 
hippocampus and cortical areas. Arch. Neurol. 57, 846-851. doi: 10.1001/arch- 
neur.57.6.846 

Holmes, A., Wrenn, C. C, Harris, A. P., Thayer, K. E., and Crawley, J. N. (2002). 
Behavioral profiles of inbred strains on novel olfactory, spatial and emo- 
tional tests for reference memory in mice. Genes Brain Behav. 1, 55-59. doi: 
10.1046/j.l601-1848.2001.00005.x 

Holsinger, R. M., Schnarr, J., Henry, P., Castelo, V. X, and Fahnestock, M. 
(2000). Quantitation of BDNF mRNA in human parietal cortex by competitive 
reverse transcription-polymerase chain reaction: decreased levels in Alzheimer's 
disease. Brain Research. Mol. Brain Res. 76, 347-354. doi: 10.1016/S0169- 
328X(00)00023-1 

Howlin, P., Goode, S., Hutton, J., and Rutter, M. (2004). Adult outcome for children 
with autism. /. Child Psychol. Psychiatry 45, 212-229. doi: 10.1 11 1/j. 1469- 
7610.2004.00215.x 

Howlin, P., and Moss, P. (2012). Adults with autism spectrum disorders. Can. J. 
Psychiatry 57, 275-283. 

Hunter, M. P., Nelson, M., Kurzer, M., Wang, X., Kryscio, R. J., Head, E., et al. 
(2011). Intersectin 1 contributes to phenotypes in vivo: implications for down 
syndrome. Neuroreport 22, 767-772. doi: 10.1097/WNR.0b013e32834ae348 

Hussain, R. J., Stumpo, D. J., Blackshear, P. J., Lenox, R. H., Abel, T, and McNamara, 
R. K. (2006). Myristoylated alanine rich C kinase substrate (MARCKS) het- 
erozygous mutant mice exhibit deficits in hippocampal mossy fiber-CA3 long- 
term potentiation. Hippocampus 16, 495-503. doi: 10.1002/hipo.20177 

Ji, L., Chauhan, A., and Chauhan, V. (2012). Reduced activity of protein kinase C 
in the frontal cortex of subjects with regressive autism: relationship with devel- 
opmental abnormalities. Int. J. Biol. Sci. 8, 1075-1084. doi: 10.7150/ijbs.4742 

Jones, E. L., Mok, K., Hanney, M., Harold, D., Sims, R., Williams, J., et al. (2013). 
Evidence that PICALM affects age at onset of Alzheimer's dementia in down 
syndrome. Neurobiol. Aging 34, 2441. el-5. doi: 10.1016/j.neurobiolaging.2013. 
03.018 

Kanner, L., Rodriguez, A., and Ashenden, B. (1972). How far can autistic children 
go in matters of social adaptation? /. Autism Child. Schizophr. 19, 213-225. 

Kelleher, R., and Bear, M. (2008). The autistic neuron: troubled translation? Cell 
135:401406. doi: 10.1016/j.cell.2008.10.017 

Khalifa, A. M., Watson-Siriboe, A., Shukry, S. G., Chiu, W. L., Nelson, M. E., Ceng, 
Y., et al. (2012). Thrl36Ile polymorphism of human vesicular monoamine 



transporter-1 (SLC18A1 gene) influences its transport activity in vitro. Neuro. 

Endocrinol. Lett. 33, 546-551. 
Liu, C. C, Kanekiyo, T., Xu, H., and Bu, G. (2013). Apolipoprotein E and Alzheimer 

disease: risk, mechanisms and therapy. Nature Rev. Neurol. 9, 106-118. doi: 

10.1038/nrneurol.2012.263 
Lohoff, F. W, Weller, A. E., Bloch, P. J., Buono, R. J., Doyle, G. A., Ferraro, T. N, 

et al. (2008). Association between polymorphisms in the vesicular monoamine 

transporter 1 gene (VMAT1/SLC18A1) on chromosome 8p and schizophrenia. 

Neuropsychobiology 57, 56-60. doi: 10.1159/000129668 
Lotter, V. (1974). Social adjustment and placement of autistic children in 

Middlesex: a follow-up study./. Autism Child. Schizophr. 4, 11-32. doi: 10.1007/ 

BF02104997 

Martin, B., Brenneman, R., Golden, E., Walent, X, Becker, K., Prabhu, V., et al. 
(2009a). Growth factor signals in neural cells: coherent patterns of interaction 
control multiple levels of molecular and phenotypic responses. /. Biol. Chem. 
284, 2493-2511. doi: 10.1074/jbc.M804545200 

Martin, B., Pearson, M., Brenneman, R., Golden, E., Wood, W. I., Prabhu, 
V., et al. (2009b). Gonadal transcriptome alterations in response to dietary 
energy intake: sensing the reproductive environment. PLoS ONE 4:e4146. doi: 
10.1371/journal.pone.0004146 

Mattson, M. P., Maudsley, S., and Martin, B. (2004). BDNF and 5-HT: a dynamic 
duo in age-related neuronal plasticity and neurodegenerative disorders. Trends 
Neurosci. 27, 589-594. doi: 10.1016/j.tins.2004.08.001 

Maudsley, S., Chadwick, W, Wang, L., Zhou, Y., Martin, B., and Park, S. S. (2011). 
Bioinformatic approaches to metabolic pathways analysis. Methods Mol. Biol. 
756, 99-130. doi: 10.1007/978-l-61779-160-4_5 

Maudsley, S., Naor, Z., Bonfil, D., Davidson, L., Karali, D., Pawson, A., et al. 
(2007). Proline-rich tyrosine kinase 2 mediates gonadotropin-releasing hor- 
mone signaling to a specific extracellularly regulated kinase-sensitive transcrip- 
tional locus in the luteinizing hormone beta-subunit gene. Mol. Endocrinol. 21, 
1216-1233. doi: 10.1210/me.2006-0053 

Mitew, S., Hay, C, Peckham, H., Xiao, J., Koenning, M., and Emery, B. (2013). 
Mechanisms regulating the development of oligodendrocytes and central ner- 
vous system myelin. Neuroscience. doi: 10. 1016/j.neuroscience.2013. 11.029. 
[Epub ahead of print] . 

Moy, S. S., and Nadler, J. J. (2008). Advances in behavioral genetics: mouse models 
of autism. Mol. Psychiatry 13, 4-26. doi: 10.1038/sj.mp.4002082 

Moy, S. S., Nadler, J. J., Young, N. B., Perez, A., Holloway, L. P., Barbaro, R. P., 
et al. (2007). Mouse behavioral tasks relevant to autism: phenotypes of 10 inbred 
strains. Behav. Brain Res. 176, 4-20. doi: 10.1016/j.bbr.2006.07.030 

Mundy, P. (2003). The neural basis of social impairments in autism: the role of 
the dorsal medial-frontal cortex and anterior cingulate system. /. Child Psychol. 
Psychiatry AA, 793-809. doi: 10.1111/1469-7610.00165 

Murray, K. D., Gall, C. M., Jones, E. G., and Isackson, P. J. (1994). Differential reg- 
ulation of brain-derived neurotrophic factor and type II calcium/calmodulin 
dependent protein kinase messenger RNA expression in Alzheimer's disease. 
Neuroscience 60, 37-48. doi: 10.1016/0306-4522(94)90202-X 

Nadler, J. J., Zou, E, Huang, H., Moy, S. S., Lauder, J., Crawley, J. N., et al. (2006). 
Large scale gene expression differences across brain regions and inbred strain 
correlate with behavioral phenotype. Genetics 174, 1229-1236. doi: 10.1534/ 
genetics.106.061481 

Nishimura, K., Lee, S. B., Park, J. H, and Park, M. H. (2012). Essential role of 
eIF5A-l and deoxyhypusine synthase in mouse embryonic development. Amino 
Acids 42, 703-710. doi: 10.1007/s00726-011-0986-z 

Noor, A., Whibley, A., Marshall, C, Gianakopoulos, P., Piton, A., Carson, A., 
et al. (2010). Disruption at the PTCHDI locus on Xp22.11 in autism spectrum 
disorder and intellectual disability. Sci. Transl. Med. 2:49ra68. doi: 10.1126/sci- 
translmed.3001267 

Park, J., Shimojo, E., and Shimojo, S. (2010). Roles of familiarity and novelty in 
visual preference judgments are segregated across object categories. Proc. Natl. 
Acad. Sci. U.S.A. 107, 14552-14555. doi: 10.1073/pnas.l004374107 

Pathania, M., Davenport, E., Muir, J., Sheehan, D., Lopez-Domenech, G., and 
Kittler, J. (2014). The autism and schizophrenia associated gene CYFIPI is criti- 
cal for the maintenance of dendritic complexity and the stabilization of mature 
spines. Transl. Psychiatry 4:e374. doi: 10.1038/tp.2014.36 

Pearson, B. L., Pobbe, R. L. H, Defensor, E. B., Osay, L., Bolivar, V. J., Blanchard, 
D. C, et al. (2011). Motor and cognitive stereotypies in the BTBR T+tf/J 
mouse model of autism. Genes Brain Behav. 10, 228-235. doi: 10.1 1 1 1601- 
183X.2010.00659.X 



Frontiers in Aging Neuroscience 



www.frontiersin.org 



September 2014 | Volume 6 | Article 225 | 24 



Jasien et al. 



Aging and autism spectrum disorder 



Petit, T. L. (1988). Synaptic plasticity and the structural basis of learning and 

memory. Neurol Neurobiol Neural. Plast. 36, 201-234. 
Phillips, H. S., Hains, J. M., Armanini, M., Laramee, G. R., and Johnson, S. A. 

(1991). BDNF mRNA is decreased in the hippocampus of individuals with 

Alzheimer's disease. Neuron 7:695702. 
Piven, J., Bailey, A. S., and Andreasen, N. (1996). Regional brain enlargement 

in autism: a magnetic resonance imaging study. /. Am. Acad. Child Adolesc. 

Psychiatry 35, 530-536. doi: 10.1097/00004583-199604000-00020 
Piven, J., and Rabins, P. (2011). Autism spectrum disorders in older adults: toward 

defining a research agenda. /. Am. Geriatr. Soc. 59, 2151-2155. doi: 10.1111/j. 

1532-5415.201 1.03632.x 
Pobbe, R. L., Defensor, E. B., Pearson, B. L., Bolivar, V. J., and Blanchard, D. 

C. (2010). Expression of social behaviors of C57BL/6J versus BTBR inbred 

mouse strains in the visible burrow system. Behav. Brain Res. 214, 443-449. doi: 

10.1016/j.bbr.2010.06.025 
Quartu, M., Lai, M. L., and Del Fiacco, M. (1999). Neurotrophin-like immunore- 

activity in the human hippocampal formation. Brain Res. Bull 48, 375-382. doi: 

10.1016/S0361-9230(99)00009-X 
Renty, J. O., and Roeyers, H. (2006). Quality of life in high-functioning 

adults with autism spectrum disorder: the predictive value of disability 

and support characteristics. Autism 10, 511-524. doi: 10.1177/1362361306 

066604 

Romanczyk, T. B., Weickert, C. S., Webster, M. J., Herman, M. M., Akil, M., 
and Kleinman, J. E. (2002). Alterations in TrkB mRNA in the human pre- 
frontal cortex throughout the lifespan. Eur. J. Neurosci. 15, 269-280. doi: 
10.1046/j.0953-816x.2001.01858.x 

Rumsey, J. M., Rapoport, J. L., and Sceery, W. R. (1985). Autistic children as adults: 
psychiatric, social, and behavioral outcomes. /. Am. Acad. Child Psychiatry 24, 
465473. doi: 10.1016/S0002-7138(09)60566-5 

Rutter, M., Greenfeld, D., and Lockyer, L. (1967). A five to fifteen year follow- 
up study of infantile psychosis (II. Social and behavioural outcome). Br. J. 
Psychiatry 113, 1183-1199. doi: 10.1 192/bjp.l 13.504. 1 183 

Scattoni, M. L., Gandhy, S. U., Ricceri, L., and Crawley, J. N. (2008). Unusual reper- 
toire of vocalizations in the BTBR T+tf/J mouse model of autism. PLoS ONE 
3:e3067. doi: 10.1371/journal.pone.0003067 

Scattoni, M. L., Ricceri, L., and Crawley, J. N. (2011). Unusual repertoire 
of vocalizations in adult BTBR T+tf/J mice during three types of social 
encounters. Genes Brain Behav. 10:4446. doi: 10.1111/j. 1601-183X.2010. 
00623.x 

Seltzer, M. M., Krauss, M. W., Shattuck, P. T, Orsmond, G., Swe, A., and Lord, C. 
(2003). The symptoms of autism spectrum disorders in adolescence and adult- 
hood. /. Autism Dev. Disord. 33, 565-581. doi: 10.1023/B:JADD.0000005995. 
02453.0b 

Shattuck, P. T, Seltzer, M. M., Greenberg, J. S., Orsmond, G. I., Bolt, D., Kring, S., 
et al. (2007). Change in autism symptoms and maladaptive behaviors in ado- 
lescents and adults with an autism spectrum disorder. /. Autism Dev. Disord. 37, 
1735-1747. doi: 10.1007/sl0803-006-0307-7 

Sokol, D. K., Maloney, B., Long, J. M., Ray, B., and Lahiri, D. K. (2011). 
Autism, Alzheimer disease, and fragile X: APP, FMRP, and mGluR5 
are molecular links. Neurology 76:13441352. doi: 10.1212/WNL.0b013e318 
2166dc7 

Stamova, B., Tian, Y., Nordahl, C, Shen, M., Rogers, S., Amaral, D., et al. 
(2013). Evidence for differential alternative splicing in blood of young 
boys with autism spectrum disorders. Mol. Autism 4:30. doi: 10.1186/2040- 
2392-4-30 



Stumpo, D. J., Bock, C. B., Tuttle, J. S., and Blackshear, P. J. (1995). MARCKS defi- 
ciency in mice leads to abnormal brain development and perinatal death. Proc. 
Natl. Acad. Sci. U.S.A. 92, 944-948. 

Suckow, A. T, Craige, B., Faundez, V., Cain, W. J., and Chessler, S. D. (2010). 
An AP-3dependent mechanism drives synaptic-like microvesicle biogenesis in 
pancreatic islet beta-cells. Am. J. Physiol Endocrinol. Metab. 299, E23-E32. doi: 
10.1152/ajpendo.00664.2009 

Szatmari, P., Bartolucci, G., Bremner, R., Bond, S., and Rich, S. (1989). A follow up 
study of high-functioning autistic children. /. Autism Dev. Disord. 19, 213-225. 
doi: 10.1007/BF02211842 

Szatmari, P., Archer, L., Fisman, S., Streiner, D. L., and Wilson, F. (1995). Asperger's 
syndrome and autism: differences in behavior, cognition, and adaptive func- 
tioning. /. Am. Acad. Child Adolesc. Psychiatry, 34, 1662-1671. doi: 10.1097/ 
00004583-199512000-00017 

Szatmari, P. (1998). Differential diagnosis of Asperger's disorder. Asperger Syndr. 
High Funct. Autism? 7, 61-76. doi: 10.1007/978-l-4615-5369-4_4 

Tapia-Arancibia, L., Aliaga, E., Silhol, M., and Arancibia, S. (2008). New insights 
into brain BDNF function in normal aging and Alzheimer disease. Brain Res. 
Rev. 59, 201-220. doi: 10.1016/j.brainresrev.2008.07.007 

Weber, A., Kohler, A., Hahn, A., Neubauer, B., and Muller, U. (2013). Benign infan- 
tile convulsions (IC) and subsequent paroxysmal kinesigenic dyskinesia (PKD) 
in a patient with 16pl 1.2 microdeletion syndrome. Neurogenetics 14, 251-253. 
doi: 10.1007/S10048-013-0376-7 

Weimer, J. M., Yokota, Y, Stanco, A., Stumpo, D. J., Blackshear, P. J., and Anton, E. 
S. (2009). MARCKS modulates radial progenitor placement, proliferation and 
organization in the developing cerebral cortex. Development 136, 2965-2975. 
doi: 10.1242/dev.036616 

Wohr, M., Roullet, F. I., and Crawley, J. N. (2011). Reduced scent marking and 
ultrasonic vocalizations in the BTBR T+tf/J mouse model of autism. Genes Brain 
Behav. 10, 35-43. doi: 10.1111/j.l601-183X.2010.00582.x 

Zhao, L., Yao, J., Mao, Z., Chen, S., Wang, Y, and Brinton, R. D. (2010). 17B- 
estradiol regulates insulin-degrading enzyme expression via an ERp7PI3-K 
pathway in hippocampus: relevance to Alzheimer's prevention. Neurobiol. Aging 
32, 1949-1963. doi: 10.1016/j.neurobiolaging.2009.12.010 

Zhubi, A, Chen, Y, Dong, E., Cook, E., Guidotti, A., and Grayson, D. (2014). 
Increased binding of MeCP2 to the GAD1 and RELN promoters may be 
mediated by an enrichment of 5-hmC in autism spectrum disorder (ASD) 
cerebellum. Transl. Psychiatry 4:e349. doi: 10.1038/tp.2013.123 

Conflict of Interest Statement: The authors declare that the research was con- 
ducted in the absence of any commercial or financial relationships that could be 
construed as a potential conflict of interest. 

Received: 27 June 2014; accepted: 08 August 2014; published online: 01 September 
2014. 

Citation: Jasien JM, Daimon CM, Wang R, Shapiro BK, Martin B and Maudsley S 
(2014) The effects of aging on the BTBR mouse model of autism spectrum disorder. 
Front. Aging Neurosci. 6:225. doi: 10.33S9/fnagi.2014.00225 
This article was submitted to the journal Frontiers in Aging Neuroscience. 
Copyright © 2014 Jasien, Daimon, Wang, Shapiro, Martin and Maudsley. This is an 
open-access article distributed under the terms of the Creative Commons Attribution 
License (CC BY). The use, distribution or reproduction in other forums is permitted, 
provided the original author(s) or licensor are credited and that the original publica- 
tion in this journal is cited, in accordance with accepted academic practice. No use, 
distribution or reproduction is permitted which does not comply with these terms. 



Frontiers in Aging Neuroscience 



www.frontiersin.org 



September 2014 | Volume 6 | Article 225 | 25 



